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14.  ABSTRACT 


This  report  results  from  a  contract  tasking  Imperial  College  London  as  follows:  1 .0  Introduction  and  Overview 


Current  fuel  quality  regulations  for  aviation  kerosene  provide  significant  scope  for  large  regional  and  source  differences.  The  chemistry  of 
high-performance  fuels  aimed  at  specific  applications  and  devices  show  less  variation  but  their  chemistries  remain  complex.  There  is  now 
arguably  a  well-established  need  to  incorporate  accurate  chemical  mechanisms  into  calculation  procedures  [1].  Flow  timescales  associated 
with  high  performance  propulsion  devices  also  increasingly  lead  to  difficulties  associated  with  the  burning  characteristics  of  the  current 
generation  of  propulsion  fuels.  In  particular,  chemical  time-scale  restrictions  pose  a  particular  challenge  in  terms  of  combustion  initiation 
and/or  flame  stability.  The  current  development  of  alternative  combustor  designs  based  premixed  or  partially  premixed  lean  burn  technologies 
further  emphasises  the  need  to  consider  finite  rate  chemistry  effects.  Critical  steps  in  the  derivation  of  realistic  chemical  mechanisms  include 
the  definition  of  a  surrogate  (or  model)  fuel  and  validation  of  the  performance  of  the  surrogate  against  practical  fuel  samples.  The  required 
practical  parameter  space  typically  leads  to  a  need  to  determine  heat  release  characteristics,  burning  velocities,  effects  of  strain,  ignition  delay 
times,  the  dynamics  of  the  fuel  breakdown  process  and  quantification  of  key  intermediate  species.  The  latter  is  often  essential  in 
understanding  the  propensity  of  a  device  to  form  pollutants,  such  as  soot,  that  may  affect  signatures  and  combustor  life.  Furthermore,  current 
indications  suggest  that  fuel  sources  will  become  significantly  more  diverse  in  the  future  and  may,  for  example,  encompass  Fischer-Tropsch 
and/or  bio-derived  components.  The  current  proposal  outlines  a  route  towards  surrogate  fuel  mechanisms  of  sufficient  accuracy  and 
generality  to  support  the  development  of  practical  devices. 

The  combination  of  advanced,  laser  based  diagnostic  techniques  and  detailed  kinetic  modelling  has  made  invaluable  contributions  to  the 
understanding  of  the  chemistry  of  fuel  conversion  [2],  Such  research  efforts  rely  upon  the  continued  development  and  application  of  innovative 
diagnostic  techniques  and  advances  in  the  understanding  of  fundamental  chemical  reaction  sequences.  Much  effort  has  been  devoted  to  the 
combustion  of  natural  gas  based  fuel  mixtures  and  past  key  differences  [3-5]  affecting  the  C1-C2  chemistry  are  now  significantly  reduced 
[6,7],  Progress  has  to  a  large  extent  been  based  on  identification  of  critical  reaction  paths  and  the  availability  of  increasingly  accurate 
determination  of  key  rate  constants  combined  with  high-quality  validation  data.  Progress  for  complex  aviation  fuel  mixtures  has  been  slower 
due  to  uncertainties  in  major  reaction  pathways,  particularly  for  aromatic  species,  limited  chemical  kinetic  and  thermochemical  data  and 
sparse  data  sets  of  limited  suitability  for  detailed  validation.  However,  the  only  approach  that  has  the  potential  to  achieve  the  objective  of 
ensuring  a  priori  knowledge  of  combustion  properties  of  a  particular  fuel  blend  is  through  the  application  of  detailed  kinetic  modelling  following 
the  independent  validation  of  the  chemistry  of  each  fuel  component  present  in  the  surrogate  mixture. 


Early  efforts  aimed  at  producing  surrogate  fuel  models  for  endothermic  and  aviation  fuels  were  pursued  by  Lindstedt  and  Maurice  [8-10]  who 
primarily  considered  n-heptane  [8],  toluene  [9],  n-decane  [10]  and  mixtures  thereof  [9,10].  The  work  resulted  a  set  of  first  generation  models 
that  identified  some  key  issues  in  the  aromatic  and  PAH  chemistry  subsequently  pursed  further  by  Lindstedt  et  al.  [11]  and  Lindstedt  and 
Rizos  [12],  The  latter  work  resulted  in  a  reassessment  of  principal  reaction  paths  leading  to  second  ring  formation.  More  recent  efforts  in  the 
European  Community  have  resulted  in  surrogate  fuel  models  based  on  two-component  mixtures  of  n-decane  and  n-propyl-benzene  [13].  The 
number  of  fuel  components  was  intentionally  limited  due  to  the  constraints  imposed  by  the  practical  application  of  such  models.  The  removal 
of  the  limitations  of  high  Damkohler  number  based  closures  for  the  chemical  reaction  rate  source  term  are  essential  in  the  context  of  a  wide 
range  of  practical  applications  where  strong  direct  kinetic  effects  are  present.  Promising  approaches  that  are  not  based  on  a  presumed  flame 
structure  include  hybrid  transported  probability  density  function/finite  volume  and  Large  Eddy  Simulation  (LES)  based  approaches.  For 
example,  Kuan  and  Lindstedt  [14]  computed  transient  turbulent  flames  with  a  large  scalar  space  featuring  20  independent  scalars  and  Sheiki 

[15]  has  highlighted  the  progress  made  in  the  application  of  LES  with  finite  rate  chemistry.  The  introduction  of  thermochemistry  into  such 
computations  requires  the  application  of  rigorous  simplification  procedures  to  reduce  the  number  of  independent  scalars. 

The  group  at  Imperial  College  [16-18]  have  further  developed  a  reaction  class  based  approach  to  model  the  oxidation  of  naphthalene,  1- 
methyl  naphthalene  and  n-propyl  benzene.  For  example,  reaction  classes  for  higher  aromatics,  such  as  naphthalene  and  indene,  were 
defined  based  on  similarities  with  the  oxidation  of  cyclopentadiene  and  benzene  [11,12],  The  application  of  the  approach  to  n-propyl  benzene 
was  shown  to  lead  to  encouraging  agreement.  The  application  of  a  surrogate  fuel  model  based  on  n-decane  and  n-propyl  benzene  to  model 
the  combustion  of  a  practical  Jet-A  sample  (collected  at  a  major  European  airport)  provided  a  further  indication  that  comparatively  simple 
surrogate  fuel  models  can  reproduce  key  oxidation  features  of  practical  fuels.  Further  information  is  available  elsewhere  [16]. 

The  work  proposed  will  serve  to  provide  consistent  reaction  mechanisms  for  surrogate  aviation  fuels  and  their  components.  To  this  effect, 
collaboration  with  USC  (Egolfopoulos,  Wang),  AFRL  (Edwards,  Williams)  and  the  wider  surrogate  fuels  working  group  and  collaborators  at 
Stanford  University  (Hanson,  Davidson). 
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2.0  Expected  Results 

1 .  A  JAVA  tool  for  the  generation  of  chemical  mechanisms  based  on  reaction  classes  relevant  to  surrogate  fuel  components. 

2.  Defined  “target”  parameters  for  the  evaluation  of  the  fidelity  of  derived  reaction  mechanisms.  Task  performed  with  project  partners. 

_3. _ Definitions  of  what  constitutes  acceptable  agreement  for  detailed  and  reduced  reaction  mechanisms  in  the  context  of  “target” 


parameters.  Task  performed  with  project  partners. 


4.  Detailed  reaction  mechanisms  derived  for  alkyl-substituted  aromatics  such  as  toluene,  m-xylene,  1 -methyl  naphthalene  and  n- 
propyl  benzene.  Revision  may  be  required  with  respect  to  selected  target  molecules  in  light  of  findings  obtained  during  the  project. 

5.  Detailed  reaction  mechanisms  derived  for  selected  n-alkane,  cyclo-paraffin  and  iso-paraffin  components. 

6.  Detailed  reaction  mechanisms  for  surrogate  fuels  evaluated  against  “target”  data.  At  least  two  fuel  components  to  be  considered  in 
the  smallest  mechanism. 


3.0  Scope  of  Activities 

All  tasks  to  be  performed  by  Imperial  College.  Collaboration  predominantly  with  AFRL,  University  of  Southern  California  (USC)  and  with  the 
Surrogate  Fuels  Working  Group.  The  tasks  outlined  below  are  sequential  and  describe  quarterly  targets  beginning  3  months  after  the 
commencement  of  the  project.  Deliverables  are  in  terms  of  brief  progress  reports,  contributions  to  EOARD/AFOSR  sponsored  meetings  and 
publications  as  appropriate.  In  addition,  reporting  will  be  performed  in  accordance  with  EOARD  and  AFOSR  requirements. 

Task  1.  Extend  existing  JAVA  tool  for  the  generation  of  reaction  mechanisms  to  incorporate  preliminary  reaction  classes  for  alkyl-substituted 
aromatics. 

Task  2.  Evaluation  of  derived  reaction  mechanisms  via  comparisons  with  experimental  data  from  the  literature  and  “target”  data  as  obtained 
by  project  partners  (USC  and  others). 

Task  3.  Extend  the  JAVA  based  tool  to  incorporate  alkane  components  including  simplified  low  temperature  chemistry. 

Task  4.  Evaluation  of  the  derived  reaction  mechanism  for  a  pure  n-alkane  component  as  well  as  for  a  preliminary  two-component  surrogate. 
Comparisons  with  agreed  “target”  parameters  and  test  cases  (AFRL,  Surrogate  Fuels  Working  Group). 

Task  5.  Refinement  of  reaction  classes  for  alkyl-substituted  aromatics  based  on  joint  work  with  project  partners  (USC). 

Task  6.  Refinement  of  reaction  classes  for  n-alkane  components  to  incorporate  a  simplified  low  temperature  mechanism  based  on 
characteristic  reaction  steps. 

Task  7.  Extend  existing  JAVA  tool  to  incorporate  preliminary  reaction  classes  for  cyclo-alkanes  derived  in  collaboration  with  project  partners 
(USC). 

Task  8.  Evaluation  of  the  derived  reaction  mechanism  for  a  selected  cyclo-alkane  via  comparisons  with  experimental  data  from  the  literature 
and  “target”  data  as  obtained  by  project  partners  (AFRL,  USC  and  others). 

Task  9.  Extend  existing  JAVA  based  tool  to  incorporate  preliminary  reaction  classes  for  iso-alkanes  derived  in  collaboration  with  project 
partners  (USC). 

Task  10.  Evaluation  of  the  derived  reaction  mechanism  for  a  selected  iso-alkane  via  comparisons  with  experimental  data  from  the  literature 
and  as  obtained  by  project  partners  (USC  and  others). 

Task  11.  Formulation  and  evaluation  of  surrogate  fuel  mechanisms  based  on  derived  and  validated  component  mechanisms.  Decision  on 
desired  dimensionality  of  the  surrogate  fuel  mechanism  (AFRL,  USC  and  the  Surrogate  Fuels  Working  Party). 

Task  12.  Final  Report 


4.0  Technical  Approach  and  Methodology 

The  current  technical  approach  is  based  on  the  further  development  of  reaction  classes  for  surrogate  fuel  components.  In  addition  to  the 
examples  mentioned  in  the  introduction,  the  technique  has  been  shown  to  produce  encouraging  results  for  aromatic  molecules  such  as 
naphthalene  and  1 -methyl  naphthalene  [16-18].  The  latter  is  a  key  surrogate  fuel  component  for  Diesel  fuel  mixtures  and  a  potential  target  fuel 
component  identified  by  the  Surrogate  Fuels  working  group.  A  further  benefit  with  the  approach  is  that  improved  chemical  kinetic  data  can 
directly  impact  related  fuel  compounds.  To  support  the  rapid  evaluation  of  mechanisms  in  light  of  improved  chemical  kinetic  and  validation 
data,  a  JAVA  based  reaction  class  tool  has  been  developed  and  applied  to  heterogeneous  reaction  systems.  The  simplification  to  encompass 
gas  phase  reactions  is  procedural.  An  important  advantage  is  that  the  approach  can  readily  be  applied  to  identify  key  sensitivities  that  merit 
detailed  computational  chemistry  and/or  experimental  investigations.  Emerging  combustion  technologies  place  greater  emphasis  on  the  low  to 
intermediate  temperature  oxidation  window  and  the  inclusion  of  appropriate  reaction  classes  has  a  tendency  to  result  in  excessive  mechanism 
sizes.  It  is  therefore  argued  that  the  approach  be  constrained  to  selected  characteristic  adducts  representative  of  wider  reaction  classes.  A 
task-by-task  outline  of  the  current  computational  project  is  given  below. 

Tasks  1,  3,  7  and  9.  A  JAVA  based  tool  has  been  created  for  the  automatic  generation  of  heterogeneous  reaction  mechanisms,  which  will  be 
modified  to  consider  homogeneous  gas  phase  reaction  classes  of  direct  relevance  to  surrogate  fuel  mechanisms.  The  initial  reaction  classes 
and  rate  data  will  be  based  on  work  performed  at  Imperial  College.  Refinements  will  take  into  account  work  performed  by  project  partners. 
The  generated  mechanisms  will  be  constrained  with  respect  to  size  by  the  selection  of  preferred  fuel  components  and  oxidation  routes. 
Further  reductions  in  size  will  result  from  path  and  sensitivity  analyses  performed  as  part  of  other  tasks. 

Tasks  2,  4,  8,  10  and  11.  Evaluation  of  derived  reaction  mechanisms  via  comparisons  with  experimental  data  from  the  literature  and  “target” 
data  as  obtained  by  project  partners.  The  reaction  mechanisms  will  also  be  subjected  to  path  and  sensitivity  analyses  with  a  view  to  reducing 
the  dimensionality  of  the  system. 

Tasks  5  and  6.  Refinement  of  reaction  classes  will  take  place  based  on  collaboration  with  project  partners  as  appropriate. 


15.  SUBJECT  TERMS 

EOARD,  Fuels,  Combustion 


16.  SECURITY  CLASSIFICATION  OF: 
a.  REPORT  I  b.  ABSTRACT  I  c. 


THIS  PAGE 


17.  LIMITATION  OF 
ABSTRACT 

UL 


18,  NUMBER 
OF  PAGES 


19a.  NAME  OF  RESPONSIBLE  PERSON 

SURYASURAMPUDI 


UNCLAS 

UNCLAS 

UNCLAS 

85 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

+44  (0)1895  616021 

Standard  Form  298  (Rev.  8/98) 
Prescribed  by  ANSI  Std.  Z39-18 


Imperial  College 

London 


Detailed  and  Simplified  Chemical  Kinetics  of  Aviation 
Fuels  and  Surrogates 

Final  Report  on  Grant  FA8655-06-1-3052 

R.P.  Lindstedt  and  V.  Markaki 

Department  of  Mechanical  Engineering,  Imperial  College  London,  Exhibition 
Road,  London,  SW7  2AZ,  UK 


Principal  Investigator  R.  P.  Lindstedt 

Submitted  to  the  European  Office  of  Aerospace  Research  and 
Development,  86  Blenheim  Crescent,  Ruislip  HA4  7HB,  UK. 

September  2009 


2 


FinalRe^orton  Grant  FA8655-06- 1 -3052 


Summary 

The  chemistries  of  aviation  fuels  are  invariably  complex  due  to  the  constituent  large 
hydrocarbon  molecules.  There  are  also  large  compositional  variations  for  a  given  fuel  type. 
Furthermore,  flow  timescales  encountered  in  high-performance  propulsion  devices  lead  to 
difficulties  associated  with  kinetically  controlled  or  influenced  phenomena  such  as  flame 
stability,  extinction  and  re-light.  Indications  are  that  fuel  sources  will  become  significantly 
more  diverse  in  the  future  and  may  encompass  Fischer-Tropsch  and/or  bio-derived 
components.  The  combustion  properties  of  such  fuels  can  vary  significantly  from  those  in 
current  use  and  this  work  outlines  a  route  towards  surrogate  fuel  mechanisms  of  sufficient 
accuracy  and  generality  to  support  the  development  of  practical  devices. 

A  reaction  class  based  route  to  the  derivation  of  detailed  chemical  kinetic  mechanisms  for 
alkyl-substituted  aromatics  is  outlined  and  applied  to  the  cyclopentadiene/indene, 
benzene/naphthalene,  toluene/ 1 -methyl  naphthalene  systems.  Work  has  also  been  extended 
to  the  n-propyl  benzene  system.  Important  reaction  channels  were  identified  for  all  systems 
and  for  selected  pathways  potential  energy  surfaces  were  determined  via  density  functional 
theory  (DFT)  and  composite  quantum  mechanical  computations.  Both  the  Rice- 
Ramsperger-Kassel-Marcus/master  equation  (RRKM/ME)  and  variable  transition  state 
theory  (VTST)  approaches  were  used  to  derive  estimates  of  the  rate  constants.  Improved 
thermodynamic  data  for  a  wide  range  of  intermediate  species  was  also  determined  at  the 
G3MP2B3  level  in  combination  with  DFT  analysis  for  internal  rotations.  The  use  of  1- 
methyl  naphthalene  as  a  surrogate  fuel  component  has  been  identified  in  the  context  of 
modulating  the  sooting  tendencies  of  aviation  surrogates  and  significant  progress  with 
respect  to  the  identification  of  oxidation  and  formation  channels  is  reported. 

Results  obtained  from  chemical  kinetic  modelling  of  cyclopentadiene,  toluene,  n-propyl 
benzene,  naphthalene  and  1 -methyl  naphthalene  oxidation  in  shock  tubes,  jet- stirred  and 
plug-flow  reactors  at  various  sets  of  representative  stoichiometries  and  temperatures 
suggest  that  the  developed  reaction  mechanims  exhibit  significantly  improved  generality. 


Acknowledgements 

The  authors  wish  to  express  their  gratitude  to  the  Air  Force  Office  of  Scientific  Research 
(AFOSR)  and  the  European  Office  of  Aerospace  Research  and  Development  (EOARD)  for 
providing  the  financial  support  for  this  work.  We  would  also  like  to  express  our  gratitude  to 
Dr.  J.  Tishkoff  and  Dr.  S.  Surampudi.  The  contributions  of  Dr.  R.  K.  Robinson,  who 
provided  much  of  the  thermochemical  and  reaction  rate  data  as  part  of  a  complementary 
effort,  are  also  gratefully  acknowledged. 


3 


Final  Report  on  Grant  FA8655-06- 1  -305 


1.0  Background 


The  increase  in  global  fuel  consumption  over  the  last  century  has  resulted  in  major 
socioeconomic  and  environmental  challenges.  Transportation  and,  to  an  increasing  extent, 
aviation,  is  one  of  the  major  fuel  consumption  sectors.  This  has  led  to  the  need  for  engines 
with  operating  envelopes  that  can  comply  with  environmental  and  energy  saving  directives. 
New  aircraft  engines  are  expected  to  perform  with  maximum  combustion  efficiency  as  well 
as  providing  stability  and  low  emissions.  Furthermore,  it  can  be  expected  that  new  fuel 
types  will  become  necessary  and  a  sound  understanding  of  current  products  forms  an 
essential  starting  point. 

Aromatic  hydrocarbons  and,  to  a  much  lesser  extent,  polycyclic  aromatic  hydrocarbons 
(PAH)  are  key  components  of  current  aviation  fuels.  It  is  well  known  that  aromatic 
compounds,  apart  from  being  responsible  for  pollutant  emissions,  also  have  carcinogenic 
and  mutagenic  properties  [1-3].  They  also  contribute  to  the  detection  of  aircraft  due  to  the 
associated  infrared  signals  produced.  The  aromatic  components  of  aviation  fuels  can  also 
reduce  the  life  cycle  of  the  combustor  as  the  associated  particulate  formation  increases  the 
radiative  heat  transfer  to  the  combustor  walls  [4]  and  also  influences  maintenance 
requirements.  Hence,  the  optimisation  of  combustor  behaviour  requires  a  good  knowledge 
of  the  chemistry  of  aromatic  fuel  components  as  well  as  calculation  methods  for  turbulent 
flows  that  are  capable  of  exploiting  such  advances. 

A  typical  aviation  fuel  consists  of  various  classes  of  hydrocarbon  compounds  and  exhibit 
different  behaviours  due  to  refining  processes  and  crude  oil  feedstocks  [5].  In  general,  the 
characteristics  of  an  aviation  fuel  is  determined  by  operational  needs.  Fuels  are  typically 
developed  to  have  good  combustion  characteristics  combined  with  good  physical  properties 
such  as  low  temperature  fluidity.  The  blends  also  need  to  consist  of  storable  hydrocarbon 
compounds.  In  the  1940’ s,  the  U.S.  Air  Force  used  ‘wide-cut’  fuel,  which  was  available  in 
large  quantities  at  that  time.  Wide-cut  fuels  were  highly  volatile  and  were  replaced  by 
kerosene  fuels  in  the  1970’s  (Jet  A,  Jet  A-l  and  JP-8)  for  safety  reasons  [6]. 

The  complexities  of  jet  fuels  are  such  that  the  direct  simulation  of  their  combustion 
behaviour  is  not  possible  due  to  the  thousands  of  chemical  compounds  present.  However, 
recent  scientific  advances  in  the  modelling  of  combustion  systems  have  provided  important 
insight  that  is  complementary  to  experimental  studies.  Such  detailed  modelling  can  be 
expected  to  provide  essential  understanding  of  the  overall  behaviour  of  a  combustion 
system  provided  that  the  chemistry  of  the  fuel  and  the  formation  of  pollutants  are 
sufficiently  well  understood.  Due  to  the  complexities  of  practical  fuels  it  is  necessary  to 
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represent  actual  fuel  mixtures  with  functionally  similar  blends  featuring  a  more  limited 
number  of  representative  compounds.  Such  surrogate  mixtures  can  be  characterized 
accurately  and  makes  it  feasible  to  study  combustion  processes  in  detail  using 
computational  methods  following  the  application  of  rigorous  systematic  reduction  methods. 

A  surrogate  aviation  fuel  can  be  chosen  to  reproduce  physical  and  chemical  properties  of  an 
aviation  fuel  such  as  heat  capacity,  enthalpy,  viscosity,  rates  of  reaction  or  specific  ignition 
and  oxidation  behaviours.  Surrogate  mixtures  can  be  used  both  computationally  and 
experimentally.  A  physical  surrogate  is  a  mixture  that  can  reproduce  physical  properties 
such  as  density,  heat  capacity,  volatility  and  a  chemical  surrogate  is  a  mixture  with  a 
chemical  class  composition  that  matches  the  real  jet  fuel  and  can  reproduce  chemical 
properties  such  as  reaction  rates,  ignition  and  sooting  behaviour.  A  surrogate  mixture  that 
has  the  same  chemical  and  physical  properties  as  the  real  fuel  is  characterized  as  a 
comprehensive  surrogate  and  can  be  expected  to  match  diverse  aspects  of  the  real  fuel 
behaviour  such  the  evaporation  process  and  the  sooting  tendency  [7]. 

The  current  work  evaluates  the  use  of  a  reaction  class  based  concept  for  the  generation  of 
chemical  mechanisms  for  surrogate  fuel  compounds.  Reaction  classes  for  higher  aromatics, 
such  as  naphthalene  and  indene  are  defined  based  on  similarities  with  the  oxidation  of 
cyclopentadiene  and  benzene  and  subsequently  applied  to  model  the  oxidation  of  toluene, 
n-propyl  benzene,  naphthalene  and  1 -methyl  naphthalene.  The  current  work  extends  past 
efforts  related  to  the  aforementioned  systems.  The  selected  compounds  have  been  identified 
as  important  in  the  context  of  a  range  of  surrogate  fuel  compositions  from  gasoline  to 
aviation  fuels.  Specifically,  methyl  groups  on  aromatic  rings  (e.g.  xylenes  and  tri-methyl 
benzenes),  or  alkyl  branches  in  general,  have  been  identified  as  important  in  the  context  of 
ignition  properties  and  1 -methyl  naphthalene  may  also  be  used  to  modulate  sooting 
tendencies  in  aviation  and  Diesel  surrogates.  These  systems  therefore  constitute  a  natural 
starting  point  for  the  evaluation  of  the  current  approach.  The  composition  of  an  aviation 
fuel  may  be  determined  via  UV  spectrophotometric  analysis  of  the  concentration  of 
individual  compounds  and  a  chromatographic  analysis  of  the  total  aromatic  content.  An 
example  of  a  chromatograph  of  a  fuel  sample  is  shown  in  Figure  1.1. 

The  construction  of  a  chemical  mechanism  is  a  complex  task  that  can  be  simplified  with  the 
identification  of  principal  reaction  pathways  and  reaction  classes  combined  with  major 
intermediate  species  that  are  formed  through  oxidation  or  pyrolysis  of  higher  hydrocarbon 
molecules.  Some  reactions  have  been  studied  experimentally  under  combustion  conditions 
and  determined  rates  can  often  be  implemented  directly  in  the  relevant  mechanism 
following  a  critical  appraisal.  Such  reactions  can  also  be  utilised  as  template  reaction  steps 
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for  the  estimation  of  relevant  routes  and  reaction  rates  of  higher  hydrocarbons  by  making 
adjustments  based  on  molecular  properties.  When  suitable  experimental  data  is  not 
available,  other  ways  of  determining  new  steps  and  reaction  rates  are  through  the 
application  of  collision  theory  [11],  quantum  mechanical  methods  [12]  and  semi-empirical 
methods  [13]. 


92 


Figure  1.1.  Chromatogram  of  a  fuel  sample,  a)  naphthalene,  b)  methyl  naphthalene,  c)  ethyl 
naphthalenes,  d)  C3-naphthalenes,  e)  diphenylamine,  f)  C4-naphthalenes,  g)  phenanthrene,  h) 
anthracene,  i)  methyl  anthracene  [9] 


Reaction  rates  that  have  been  determined  with  all  of  the  aforementioned  methods  are 
utilised  in  the  present  work  and  details  are  provided  for  each  specific  system  in  the 
following  sections.  Important  reaction  channels  were  identified  for  all  systems  and  for 
selected  pathways  potential  energy  surfaces  (PES)  were  determined  via  density  functional 
theory  (DFT)  and  composite  quantum  mechanical  computations  [15].  Both  the  Rice- 
Ramsperger-Kassel-Marcus/master  equation  (RRKM/ME)  and  variable  transition  state 
theory  (VTST)  approaches  were  used  to  derive  estimates  of  the  rate  constants.  Improved 
thermodynamic  data  for  a  wide  range  of  intermediate  species  was  also  determined  at  the 
G3MP2B3  level  in  combination  with  DFT  analysis  for  internal  rotations  [15].  Data  has  also 
been  obtained  from  Burcat  and  Ruscic  [14]  where  appropriate.  The  parts  of  the  detailed 
chemical  mechanisms  discussed  in  the  current  report  can  be  found  in  Appendix  A. 
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2.0  The  Cyclopentadiene  System 

The  study  of  the  chemical  kinetics  of  aromatic  hydrocarbons  in  combustion  processes  has 
become  a  main  focus  of  combustion  research.  Apart  from  the  environmental  impact  caused 
by  soot  formation,  polyaromatic  compounds  have  a  biological  impact  on  human  health. 
However,  detailed  chemical  models  are  often  tentative.  There  is  a  limited  amount  of 
elementary  reaction  rate  data  in  the  literature,  which  is  mostly  focussed  on  single  ring 
molecules  such  as  benzene  and  toluene.  Benzene  oxidizes  to  the  phenoxy  radical  which 
decomposes  to  the  cyclopentadienyl  radical  and  carbon  monoxide.  The  cyclopentadienyl 
radical  is  subject  to  ring  opening  and  forms  a  transition  point  between  aromatic  and 
aliphatic  compounds  [24].  Cyclopentadiene  is  also  an  important  intermediate  species  in  the 
combustion  of  other  single  ring  aromatics.  The  potential  importance  of  the  cyclopentadienyl 
radical  in  PAH  growth  has  also  been  noted  in  various  studies  [25-31].  Hence,  an  adequate 
knowledge  of  the  cyclopentadienyl  combustion  chemistry  is  essential  for  the  accurate 
modelling  of  the  aromatic  components  of  the  aviation  fuels. 

The  pyrolysis  of  cyclopentadiene  has  been  extensively  studied  behind  reflected  shock 
waves  [24,  32-36],  but  there  is  only  one  detailed  oxidation  study,  performed  by  Butler  and 
Glassman  [24]  using  a  plug  flow  reactor.  It  is  generally  agreed  that  the  cyclopentadiene 
consumption  under  pyrolytic  conditions  begins  with  CH  fission  leading  to  the  formation  of 
the  cyclopentadienyl  radical  plus  a  hydrogen  atom  (13). 

C5H6  =  C5H5  +  H  (13) 

There  are  significant  differences  in  reaction  rates  obtained  in  the  different  studies.  Roy  et 
al  [34]  studied  the  pyrolysis  of  cyclopentadiene  in  argon  mixtures  over  a  pressure  range  of 
0.7  to  5.6  bar  and  obtained  a  rate  constant  of  k13=  4.0  x  1014  exp  (-38760/T)  s"1  for  which 
they  concluded  that  there  was  no  pressure  dependence.  Kern  et  al  [36]  studied 
cyclopentadiene  pyrolysis  at  reduced  pressures  in  the  range  100-450  Torr  and  proposed  a 
barrier  of  368  kJ/mol  for  the  CH  fission  when  combined  with  a  non-Arrhenius  fit.  Colket  et 
al.  [35]  proposed  a  rate  constant  of  k13=  2  x  1015  exp(-339  kJ  mol'VRT)  s"1  for  a  range  of 
slightly  higher  pressures  than  the  study  of  Roy  et  al  [34].  Zhong  and  Bozzelli  [37] 
suggested  a  rate  of  k13=  5.96  x  1014  exp  (-314  kJ  mol'VRT)  s'1,  which  was  calculated  using 
the  Quantum  Rice-Ramsperger  Kassel  (QRRK)  approach  at  atmospheric  pressure  and  a 
temperature  range  of  900  -  1300  K. 

Roy  et  al  [33]  studied  the  reverse  reaction,  which  involves  the  cyclopentadienyl 
recombination  with  a  hydrogen  atom  (-13),  behind  reflected  shock  waves  and  a  rate 
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constant  of  k_13=2.6  x  1011  m3  kmol"1  s'1  was  obtained  at  a  pressure  of  2  bar.  Hence,  in 
combination  with  the  previous  study  by  the  same  group,  an  equilibrium  constant  for 
reaction  (13)  was  calculated.  Moreover,  data  for  the  enthalpy  of  formation  of  the  C5H5  were 
obtained.  Zhong  and  Bozzelli  [37]  also  proposed  a  high-pressure  limit  rate  constant  of  k_13= 
3.2  x  1011  m3  kmol1  s'1  for  reaction  (-13).  The  rate  proposed  was  tested  against  the  data  of 
Lovell  et  al.  [38]  at  1  atm  and  1100  K  and  there  were  discrepancies  of  about  a  factor  of  2 
between  the  model  and  measurements.  In  the  current  work  the  rate  determined  by  Kern  et 
al.  [36]  has  been  adopted. 

The  H  abstraction  leading  to  the  formation  of  the  cyclopentadienyl  radical  and  molecular 
hydrogen  is  another  important  reaction  channel  for  the  decomposition  of  cyclopentadiene. 

C5H6  +  H  =  C5H5  +  H2  (9) 

Roy  et  al.  [34]  studied  reaction  (9)  and  obtained  a  rate  constant  of  k9  =  2.8  x  1010  exp(- 
1 137/T)  m3  kmol1  s'1.  On  the  other  hand  Emdee  et  al.  [39]  proposed  a  rate  constant  of  k9  = 
2.19  x  105  x  T177  exp(-125  kJ  mol'VRT)  m3  kmol'1  s'1.  Moskaleva  and  Lin  [40]  studied  the 
reaction  between  hydrogen  and  cyclopentadiene  at  the  modified  Gaussian-2  level  of  theory 
and  provided  potential  energy  surfaces  and  reaction  rates  for  the  channel.  The  transition 
state  (TST)  and  Rice-Ramsperger-Kassel-Marcus  (RRKM)  calculations  performed 
produced  a  rate  constant  k9  =  3.03  x  105  T1 71  exp(-117  kJ  mol'VT)  m3  kmol"1  s"1  valid  for  a 
temperature  range  1000  K  <  T  <  3000  K.  The  study  showed  that  the  abstraction  rate 
constant  does  not  have  any  pressure  dependency,  but  a  positive  temperature  dependence. 
Moskaleva  and  Lin  [40]  also  showed  that,  apart  from  hydrogen  abstraction,  hydrogen 
addition  to  the  cyclopentadiene  ring  should  also  be  included  in  the  modelling  of  the  system. 

The  CH  fission  leading  to  the  cyclopentadienyl  radical  is  followed  by  C-C  breakage  and  the 
formation  of  acetylene  and  the  propargyl  radical  (3).  Studies  have  shown  that  this  reaction 
is  a  multi-step  process  which  leads  to  a  stable  intermediate  and  then  a  ring  opening 
producing  open  chain  radicals  followed  by  the  C-C  fission  for  the  formation  of  the  final 
products  [41]. 

C5H5  =  C3H3  +  C2H2  (3) 

Roy  et  al.  [34]  performed  computations  for  reaction  channel  (3)  and  optimized  the 
geometries  of  the  cyclopentadienyl  radical  and  of  all  the  intermediates,  the  transition  states 
and  products  at  the  UHF/6-31G*  level.  According  to  their  study,  the  cyclic  C5H5  decay  is 
initiated  by  a  1,2-H  transfer  followed  by  a  ring  opening  and  isomerization  to  a  number  of 
straight  chain  radicals  that  lead  to  the  formation  of  acetylene  and  the  propargyl  radical. 
However,  application  of  the  proposed  rate  resulted  in  an  over-prediction  of  acetylene 
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concentrations  by  a  factor  of  three.  Hence,  they  reduced  the  rate  constant  by  a 
corresponding  factor  and  achieved  a  good  agreement  with  measurements.  In  addition  to  this 
study,  Moskaleva  and  Lin  [40]  estimated  the  barrier  heights  of  the  1,2  hydrogen  transfer  in 
the  cyclopentadienyl  radical  (274.9  kJ/mol),  ring  opening  (317.1  kJ/mol)  and  C-C  breakage 
(317.1  kJ/mol),  which  are  in  agreement  with  values  proposed  by  Roy  et  al.  [34]  obtained 
using  lower  level  computations  [41].  Kern  et  al.  [36]  also  obtained  a  rate  constant  for 
reaction  (3)  using  RRKM  computations.  The  1,2-H  migration  was  identified  as  the  rate- 
limiting  step  with  a  barrier  of  259.0  kJ/mol.  The  result  is  in  contrast  to  Roy  et  al.  [34]  who 
proposed  the  ring  opening  as  the  rate  limiting  step.  Kern  et  al.  [36]  also  showed  that  there 
is  a  10-fold  reaction-path-degeneracy  that  arises  from  a  facile  pseudorotation  in  the 
molecule. 

Cyclopentadiene,  apart  from  the  unimolecular  decomposition  leading  to  the  formation  of 
acetylene  and  the  propargyl  radical  which  contribute  to  PAH  growth,  can  also  recombine 
with  itself  or  its  radical  leading  to  the  formation  of  naphthalene  or  indene  [26,  42].  There 
are  a  number  of  studies  that  have  exploited  this  potential  formation  mechanism.  Melius  et 
al.  [26]  studied  the  recombination  of  two  cyclopentadienyl  radicals  that  lead  to  the 
formation  of  naphthalene  via  a  nine  step  mechanism  that  goes  through  the  formation  of 
hydrofulvalene  and  involves  a  three-membered  ring  closing  and  opening  of  resonance- 
stabilized  radicals.  Kislov  and  Mebel  [43]  showed  that  at  low  temperatures  naphthalene  was 
the  major  product  whereas  at  high  temperatures  fulvalene  was  dominant.  Lu  et  al.  [44] 
extended  the  mechanism  of  PAH  growth  of  naphthalene  formation  from  cyclopentadienyl 
recombination,  as  proposed  by  Melius  et  al.  [26],  to  compounds  that  contain  four  six- 
membered  rings. 

Lindstedt  et  al.  [45]  and  Lindstedt  and  Rizos  [46]  questioned  the  global  reaction  rate 
proposed  by  Marinov  et  al.  [42],  as  computations  [46]  showed  that  naphthalene  levels  were 
up  35  times  higher  than  measurements  in  systems  with  high  cyclopentadienyl 
concentrations.  Moreover,  McEnally  and  Pfefferle  [47]  performed  experimental  work  with 
13C-labeled  aromatic  compounds  and  showed  that  the  labelled  cyclopentadienyl  moieties 
did  not  contribute  significantly  to  second  ring  formation  under  their  conditions.  Hence,  an 
alternative  two-step  reaction  featuring  stabilization  of  C5H5-C5H4  with  a  a  barrier  of  34 
kJ/mol  was  proposed  by  Lindstedt  et  al.  [45].  The  applied  barrier  thus  follows  the  potential 
energy  surface  determined  by  Melius  et  al.  [26]. 

The  potential  role  of  cyclopentadienyl  moieties  in  PAH  growth  was  also  studied  by 
Mulholland  and  co-workers  [48]  and  the  results  suggest  that  indene,  naphthalene  and 
benzene  are  major  products  of  cyclopentadiene  pyrolysis.  Apart  from  the  cyclopentadienyl 
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radical  route  to  naphthalene,  recombination  has  also  been  proposed  as  route  to  other 
aromatics.  The  importance  of  such  reaction  for  PAH  growth  was  studied  by  Violi  and  co¬ 
workers  [49-51],  who  proposed  a  two-step  molecule-radical  addition  mechanism  followed 
by  rearrangement.  Wang  et  al.  [25]  applied  the  molecule  radical  mechanism  in  order  to 
produce  new  pathways  for  the  formation  of  indene  and  naphthalene  during 
cyclopentadienyl  pyrolysis.  Four  new  reaction  pathways  for  aromatic  growth  from 
cyclopentadiene  pyrolysis  were  proposed  leading  to  naphthalene,  indene  and  benzene  via 
intramolecular  addition,  C-H  and  C-C  scissions.  Density  functional  theory  calculations 
were  performed  to  calculate  transition  states,  energy  barriers  for  isomerizations,  (3-scissions 
and  dissociation  reactions.  The  study  suggested  that  indene  was  mainly  formed  by 
intramolecular  addition  of  cyclopentadiene  to  cyclopentadienyl  via  a  resonantly  stabilised 
cyclopentadiene-cyclopentadienyl  dimer. 

Kislov  and  Mebel  [52]  also  studied  indene  formation  via  the  combination  of 
cyclopentadiene  and  the  cyclopentadienyl  radical  and  showed  that  at  temperatures  relevant 
to  combustion,  indene  was  the  major  product  (>  50%).  The  mechanism  suggested, 
combined  with  their  computed  tests,  was  in  agreement  with  experimental  data  for 
cyclopentadienyl  pyrolysis  that  showed  both  naphthalene  and  indene  as  major  products. 
Another  reaction  route  that  involves  the  cyclopentadienyl  radical  features  its  combination 
with  the  methyl  radical  leading  to  the  formation  of  fulvene  (C5H4CH2),  which  subsequently 
isomerizes  to  benzene.  Moskaleva  et  al.  [53]  studied  this  reaction  route  and  showed  that  the 
combination  of  the  two  radicals  produce  an  intermediate  specie  (C5H5CH3)  that  proceeds  to 
the  formation  of  C5H4CH3  +  H.  The  hydrogen  elimination  occurs  from  the  ring  since  the  C- 
H  bond  is  weaker  than  that  of  the  methyl  group.  The  second  hydrogen  atom  elimination 
comes  from  C5H4CH3  leading  to  the  formation  of  fulvene.  The  reaction  channel  was  also 
studied  by  Melius  et  al.  [26]  who  underlined  the  importance  of  hydrogen  migration  around 
the  cyclopentadienyl  moiety  in  providing  resonance-stabilized  radicals.  The  findings  are 
also  supported  by  a  more  recent  study  of  Lindstedt  and  Rizos  [46]  on  the  oxidation  of 
cyclopentene. 

2.1  The  Current  Contribution 

The  starting  point  for  the  current  work  stems  from  previous  studies  related  to  the  chemistry 
of  aromatics  [45,  46],  the  oxidation  of  fine  carbon  based  particles  [54]  and  an  earlier  study 
by  Zhong  and  Bozzelli  [37].  The  current  report  is  focussed  on  the  oxidation  of 
cyclopentadiene,  but  pyrolysis  cases  have  also  been  studied  with  similar  agreement  to  that 
shown  below. 
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The  developed  mechanism  was  validated  under  the  plug  flow  reactor  conditions  used  by 
Butler  [55].  Rates  were  analysed  for  both  fuel  lean  and  rich  conditions  and  for  temperatures 
varying  from  1100  to  1200  K.  The  studies  identified  possible  reaction  channels  featuring  O, 
OH,  H02  and  02,  though  estimates  of  the  rates  of  reaction  have  proved  problematic  in  the 
past.  In  the  current  work,  reaction  rates  for  isomerization  reactions  and  the  thermal 
dissociation  of  oxygenated  C5  species  were  updated  using  potential  energy  surfaces  (PES) 
determined  via  density  functional  theory  (DFT)  and  composite  quantum  mechanical 
computations  [15].  Both  the  Rice-Ramsperger-Kassel-Marcus/master  equation  approach 
(RRKM/ME)  and  variable  transition  state  theory  (VTST)  approaches  were  used  to  derive 
estimates  of  the  rate  constants  [15].  Improved  thermodynamic  data  for  a  wide  range  of 
intermediate  species  was  also  determined  at  the  G3MP2B3  level  in  combination  with 
density  functional  theory  analysis  for  internal  rotations. 

The  oxidation  of  cyclopentadiene  was  studied  under  plug  flow  reactor  conditions  obtained 
from  Butler  [55].  The  full  reaction  mechanism  applied  here  consists  of  1431  reversible 
reactions  and  269  chemical  species.  Three  oxidation  cases  were  computed  for 
stoichiometric  and  rich  fuel  mixtures  at  atmospheric  pressure  utilizing  nitrogen  as  the 
carrier  gas.  The  conditions  are  presented  in  Table  2.1. 


Case 

0 

To 

(K) 

Initial  Fuel 
(ppm) 

Initial  Oxygen 
(ppm) 

Modelling 
Time  Shift 

(ms) 

Experimental 
Time  Shift 

(ms) 

1 

1.03 

1198 

2243 

14128 

20 

36 

2 

1.03 

1148 

1051 

6618 

50 

123 

3 

1.61 

1153 

2070 

8363 

20 

50 

Table  2.1  Experimental  and  modelling  conditions  for  cyclopentadiene  oxidation  in  a  flow  reactor. 
(The  experimental  time  shifts  correspond  to  the  time  shift  applied  in  the  computations  performed  by 
Butler  [55]). 


A  reaction  rate  analysis  was  performed  for  a  stoichiometric  fuel  mixture  with  O  =  1.03  at  T 
=  1148  K,  P  =  1  atm  and  an  initial  fuel  concentration  of  1051  ppm.  The  fuel  is  consumed 
via  four  major  channels  that  involve  hydrogen  abstraction  via  H,  O  and  OH  radicals  and 
thermal  decomposition. 

Major  species  concentrations  are  compared  to  measurements  as  shown  in  Figs.  2.1  -  2.6.  It 
can  be  seen  that  there  is  reasonably  good  agreement.  However,  surface  chemistry 
interference  with  the  gas  phase  chemistry  was  identified  as  a  problem  by  Butler  [55]  and 
some  caution  is  hence  required.  Moreover,  during  the  experimental  studies,  measurements 
showed  an  initial  discontinuous  drop  in  the  fuel  concentration  at  the  first  data  point  due  to 
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the  recirculation  zone  in  the  diffuser  and  surface  chemistry  near  the  throat  of  the  reactor 
[55].  Hence,  a  time  shift  is  necessary  in  order  to  reproduce  the  species  profiles. 

Reaction  (9)  is  responsible  for  35%  of  the  total  fuel  consumption  and  is  assigned  a  rate 
proposed  by  Robinson  [15].  The  channel  is  of  major  importance  as  small  perturbations  to 
the  reaction  rate  causes  big  differences  in  the  fuel  consumption  profile.  Approximately  25% 
of  the  fuel  decay  proceeds  via  reaction  (11)  that  was  assigned  a  rate  proposed  by  Lindstedt 
and  Rizos  [46].  A  rate  discussed  by  Leung  and  Lindstedt  [58]  was  used  for  reaction  (10), 
which  consumes  a  further  10%.  Reaction  (13)  is  responsible  for  12%  of  the  total  fuel 
consumption  and  was  assigned  a  rate  proposed  by  Kern  et  al.  [36]. 


Figure  2.1  Concentration  profiles  during  cyclopentadiene  oxidation  in  a  plug  flow  reactor  for  ®  = 
1.03,  P  =1  atm,  T  =  1198  K,  initial  fuel  concentration  2243  ppm  (Case  1,  Table  2.1).  Circles  are 
measurements  [55]  and  the  solid  line  the  current  simulation. 


Figure  2.2  Concentration  profiles  during  cyclopentadiene  oxidation  in  a  plug  flow  reactor  for  ®  = 
1.03,  P  =1  atm,  T  =  1198  K,  initial  fuel  concentration  2243  ppm  (Case  1,  Table  2.1).  Circles  are 
measurements  [55]  and  the  solid  line  the  current  simulation. 
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Figure  2.3  Concentration  profiles  during  cyclopentadiene  oxidation  in  a  plug  flow  reactor  for  ®  = 
1.03,  P  =1  atm,  T  =  1148  K,  initial  fuel  concentration  1050  ppm  (Case  2,  Table  2.1).  Circles  are 
measurements  [55]  and  the  solid  line  the  current  simulation. 


Figure  2.4  Concentration  profiles  during  cyclopentadiene  oxidation  in  a  plug  flow  reactor  for  ®  = 
1.03,  P  =1  atm,  T  =  1148  K,  initial  fuel  concentration  1050  ppm  (Case  2,  Table  2.1).  Circles  are 
measurements  [55]  and  the  solid  line  the  current  simulation. 

Around  85%  of  the  fuel  consumption  leads  to  the  cyclopentadienyl  radical.  The 
consumption  of  C5H5  follows  two  major  routes.  The  pathway  leading  to  Ci0H9F  (6)  is 
responsible  for  33%  of  the  consumption  and  was  assigned  a  rate  proposed  by  Lindstedt  et 
al.  [45].  The  second  major  consumption  channel  (33%)  occurs  via  thermal  decomposition 
of  the  cyclopentadienyl  radical  via  C-C  scission  leading  to  acetylene  and  the  propargyl 
radical  (3).  The  channel  is  of  significant  importance  and  affects  the  temporal  evolution  of 
the  acetylene  concentration  profile.  A  rate  of  Kern  et  al.  [36]  was  also  evaluated,  but  found 
to  lead  to  an  overproduction  of  acetylene  by  50%.  Hence,  an  adjustment  by  a  factor  of  2 
was  applied  leading  to  more  reasonable  acetylene  profiles. 
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Figure  2.5  Concentration  profiles  during  cyclopentadiene  oxidation  in  a  plug  flow  reactor  for  ®  = 
1.61,  P  =1  atm,  T  =  1153  K,  initial  fuel  concentration  2070  ppm  (Case  3,  Table  2.1).  Circles  are 
measurements  [55]  and  the  solid  line  the  current  simulation. 


Figure  2.6  Concentration  profiles  during  cyclopentadiene  oxidation  in  a  plug  flow  reactor  for  ®  = 
1.61,  P  =1  atm,  T  =  1153  K,  initial  fuel  concentration  2070  ppm  (Case  3,  Table  2.1).  Circles  are 
measurements  [55]  and  the  solid  line  the  current  simulation. 


The  oxidation  channels  featuring  H02  attack  leading  to  the  formation  of  C5H4O  (7)  and 
C5H5O  (5)  are  each  responsible  for  8%  of  the  C5H5  consumption. 


c5H5 

+  C5H5 

=  c10H9F 

+ 

H 

(6) 

c5h5 

=  C3H3 

+ 

c2H2 

(3) 

c5h5 

+  ho2 

=  C5H40 

+ 

h2o 

(7) 

c5h5 

+  ho2 

=  C5H50 

+ 

OH 

(5) 

The  C10H9F  moiety  is  a  precursor  to  naphthalene,  which  is  subsequently  formed  via  thermal 
dissociation  (79).  Naphthalene  decomposes  through  two  major  reaction  channels  that 
involve  hydrogen  abstraction  (78)  via  OH  (49%)  and  oxygen  addition  (77)  forming  Ci0H7O 
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(39%)  that  were  assigned  rates  adopted  from  the  kinetics  of  benzene  and  adjusted  according 
to  molecular  weight  differences. 


c10h9f 

II 

n 

0 

a 

00 

+ 

H 

(79) 

CioH8 

+  OH 

=  c10H7 

+ 

H20 

(78) 

CioH8 

+  0 

=  c10Hvo 

+ 

H 

(77) 

Indene  is  formed  predominantly  (48%)  through  acetylene  recombination  with  the  benzyl 
radical  (34).  The  rate  assigned  to  this  channel  was  adopted  from  Colket  et  al.  [59].  An 
additional  44%  of  the  indene  formation  comes  from  the  recombination  of  cyclopentadiene 
and  cyclopentadienyl  radical  (12)  with  the  simultaneous  abstraction  of  a  methyl  radical. 
The  reaction  channel  was  studied  extensively  by  Wang  et  al.  [25],  who  proposed  possible 
intermediate  routes  for  the  naphthalene  and  indene  formation  via  C5H6  (CPD)  and  C5H5 
(CPDyl)  recombination.  In  this  study  it  is  assumed  that  the  recombination  of  the  two 
species  does  not  stabilize  to  the  CPD-CPDyl  intermediate,  but  leads  to  a  bridged 
intermediate  specie  with  the  radical  on  the  bridged  atom.  The  highest  energy  barrier  (177 
kJ/mol)  that  occurs  from  the  conversion  of  this  bridged  specie  to  the  final  indene  molecule 
through  other  intermediate  routes,  as  discussed  by  Wang  et  al.  [25],  was  used  as  the 
activation  energy  of  the  global  step. 

C7H7  +C2H2  =  C9H8  +H  (34) 

C5H5  +  C5H6  =  C9H8  +  CH3  (12) 

Indene  decomposes  solely  to  the  indenyl  radical  via  three  major  reaction  routes.  Reaction 
(73)  is  responsible  for  59%  of  the  indene  consumption,  (65)  for  21%  and  (74)  12% 
respectively.  Reactions  rates  utilized  for  these  three  indenyl  consumption  routes  are  adopted 
from  suggestions  by  Potter  [60]  based  on  the  kinetics  of  cyclopentadiene  with  reaction  rates 
proposed  by  Lindstedt  and  Rizos  [46]  and  Rizos  [57]. 


C9H8 

+  H 

=  c9H7 

+  h2 

(73) 

c9h7 

+  H 

=  c9H8 

(65) 

c9h8 

+  OH 

=  c9H7 

+  h2o 

(74) 

Approximately  63%  of  the  total  indenyl  radical  comes  from  indene  and  23%  is  produced 
via  CO  expulsion  from  Ci0H7O,  which  is  one  of  the  major  products  of  naphthalene 
oxidation.  The  indenyl  radical  is  the  precursor  of  phenylacetylene  via  the  formation  of 
C9H70.  The  indenyl  radical  oxidizes  to  C9H70  (83%)  via  reaction  with  the  H02  radical.  The 
reaction  rate  utilized  for  this  channel  was  obtained  by  Lindstedt  et  al.  [54].  The  C9H70 
molecule  decomposes  to  C8H7  (100%)  via  CO  expulsion  followed  by  thermal  dissociation 
of  C8H7  leading  to  the  formation  of  phenylacetylene. 
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c9h7  +  ho2 

=  c9h7o 

+  OH 

(68) 

c9h7o 

r- 

DC 

00 

O 

II 

+  CO 

(75) 

c8h7 

=  c8H6 

+  H 

(58) 

Phenylacetylene  is  consumed  via  four  major  reaction  pathways  that  predominantly  feature 
oxygen  attack.  The  displacement  of  the  chain  via  oxygen  atom  attack  is  responsible  for  33% 
of  the  phenylacetylene  consumption  (54)  and  the  step  features  a  rate  suggested  by  Lindstedt 
et  al.  [45].  The  second  major  channel  (22%)  involves  hydrogen  abstraction  via  OH  (52).  A 
rate  adopted  by  Frenklach  et  al.  [61]  was  assigned  to  this  channel.  Moreover,  reactions  (53) 
and  (51)  involve  hydrogen  abstraction/oxygen  addition  and  consume  12%  and  11%, 
respectively.  Rates  applied  to  these  channels  were  adopted  from  Potter  [60]. 


c8h6 

+  0 

=  C6H5 

+ 

c2HO 

(54) 

c8h6 

+  OH 

=  c8H5 

+ 

h2o 

(52) 

c8h6 

+  0 

=  c8h50 

+ 

H 

(53) 

c8h6 

+  0 

=  c6h5c2o 

+ 

H 

(51) 

Toluene  is  formed  from  the  benzyl  radical,  with  the  formation  of  the  benzyl  radical 
controlled  (84%)  by  the  recombination  of  the  cyclopentadienyl  radical  with  acetylene.  The 
rate  for  the  channel  originally  stems  from  the  work  of  Colket  and  Seery  [62]. 

C5H5  +  C2H2  =  C7H7  (2) 

C7H8  =  C7H7  +  H  (38) 

Approximately  22%  of  toluene  leads  to  the  formation  of  benzene  through  a  chain 
displacement  reaction  via  hydrogen.  However,  benzene  is  predominantly  formed  (94%) 
through  isomerisation  from  fulvene  (26).  The  isomerization  of  fulvene  to  benzene  is 
assigned  a  rate  proposed  by  Marinov  et  al.  [63].  Fulvene  is  mostly  produced  via 
isomerization  reactions  from  C6H6(S)  (43%)  and  C6H6(B)  (26%)  and  via  thermal  hydrogen 
dissociation  (26%)  from  C5H4CH3.  The  methyl  cyclopentadienyl  radical  is  formed  via 
recombination  of  CH3  with  C5H5  (93%)  and  was  assigned  a  rate  adopted  from  suggestions 
of  Lindstedt  et  al  [45].  Acetylene  is  predominantly  formed  (77%)  by  the  thermal 
decomposition  of  the  cyclopentadienyl  radical. 

C5H5  =  C2H2  +C3H3  (3) 

A  reaction  rate  analysis  was  also  performed  for  a  fuel  rich  mixture  with  ®  =  1.61  at  T  = 
1153  K,  P  =  1  atm  and  an  initial  fuel  concentration  of  2070  ppm.  The  fuel  consumption 
follows  the  same  routes  as  in  the  stoichiometric  case.  However,  the  impact  of  reaction  (9) 
increases  by  4%  and  the  impact  of  reactions  (11)  and  (14)  reduce  by  4-5%  compared  to  the 
stoichiometric  case. 


c5H6 

+  H 

=  c5H5 

+  h2 

(9) 

c5h6 

+  OH 

=  c5H5 

+  h2o 

(11) 

c5h7 

=  c5H6 

+  H 

(14) 

The  cyclopentadienyl  radical  follows  two  major  consumption  routes  under  fuel  rich 
conditions,  as  was  the  case  for  the  stoichiometric  conditions,  and  38%  of  C5H5  leads  to 
C10H9F  via  recombination  and  26%  thermally  decomposes  to  acetylene  and  the  propargyl 
radical.  The  reactions  featuring  H02  attack  and  leading  to  formation  of  C5H5O  and  C5H4O 
perform  similarly  to  the  stoichiometric  case  and  found  to  be  responsible  for  8%  of  total 
cyclopentadienyl  consumption. 

The  naphthalene,  indene,  styrene,  toluene,  butadiene  and  acetylene  formation  and 
consumption  routes  follow  the  same  behaviour  as  in  stoichiometric  mixtures  with  a 
perturbation  of  order  4-6%  for  each  of  the  steps.  Apart  from  the  rate  analysis  performed  for 
a  temperature  range  of  1149  -  1153  K,  a  test  was  also  performed  for  a  stoichiometric  fuel 
mixture  (Case  1  -  Table  2.1)  at  the  temperature  of  T  =  1 198  K  in  order  to  identify  important 
fuel  decomposition  pathways.  Reactions  (9)  and  (11)  were  found  to  be  dominant  and 
contribute  33%  and  27%  to  the  fuel  consumption.  Moreover,  reaction  (10)  is  responsible 
for  12%  at  the  higher  temperature. 


c5h6 

+  H 

=  c5h5 

+  h2 

(9) 

c5h6 

+  OH 

=  c5h5 

+  h2o 

(11) 

c5h6 

+  0 

=  c5h5 

+  OH 

(10) 

It  should  be  highlighted  that  once  the  cyclopentadienyl  radical  is  formed,  its  consumption  is 
controlled  by  two  dominant  reaction  channels  that  are  identical  to  the  C5H5  consumption 
route  for  rich  mixtures,  with  the  difference  that  the  role  of  reactions  (6)  and  (3)  is  somewhat 
affected.  The  impact  of  reaction  (3)  increases  from  32%  for  stoichiometric  and  26%  for  rich 
mixtures  at  T  =  1 149  K  to  34%  for  stoichiometric  mixtures  at  a  temperature  of  T  =  1 198  K. 
Whereas  the  contribution  of  the  cyclopentadienyl  recombination  reaction  falls  from  38%  for 
rich  and  33%  for  stoichiometric  mixtures  at  the  lower  temperature,  to  24%  for 
stoichiometric  mixtures  at  the  higher  temperature.  Accordingly,  the  more  the  temperature 
increases,  the  more  the  thermal  decomposition  of  C5H5  dominates  consumption. 

C5H5  +  C5H5  =  C10H9F  +  H  (6) 

C5H5  =  C3H3  +  C2H2  (3) 

Due  to  the  impact  of  the  C5H5  decomposition  channels  on  the  evolution  of  the 
cyclopentadienyl  radical,  the  concentration  of  the  naphthalene  is  increased  by  about  50% 
compared  to  the  lower  temperature  case  of  about  T  =  1 150  K.  Added  to  this,  the  fuel  at  T  = 
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1198  K  shows  a  steep  decomposition  profile  and  the  acetylene  concentration  is  double 
compared  to  the  T  =  1 153  K  oxidation  case. 

The  increased  impact  of  the  thermal  decomposition  route  of  the  C5H5,  is  responsible  for  the 
increased  propargyl  radical  and  acetylene  concentrations.  The  propargyl  radical  pool 
follows  a  recombination  route  that  produces  benzene.  Hence,  as  expected,  the  benzene 
concentration  is  increased  by  60%  at  the  temperature  of  T  =  1198  K  as  compared  to  the  T  = 
1150  K  cases.  Moreover,  the  concentration  of  benzene  is  affected  by  the  increased 
contribution  of  reaction  (1)  leading  to  C5H4CH3,  which  is  the  precursor  to  fulvene. 

C5H5  +CH3  =  C5H4CH3  +  H  (1) 

The  formation  and  consumption  routes  of  the  rest  of  the  intermediate  species,  such  as 
indene,  phenylacetylene,  toluene  follow  the  same  behaviour  as  at  lower  temperatures. 


3.0  The  Oxidation  of  Toluene 

Alkylated  benzenes  are  a  significant  hydrocarbon  class  in  practical  fuels,  hence  a  detailed 
knowledge  of  their  combustion  chemistry  is  necessary.  Toluene  is  the  simplest  example  and 
its  chemistry  constitutes  the  starting  point  for  the  development  of  mechanisms  describing 
other  alkylated  benzenes.  Specifically,  in  the  current  work  the  molecule  forms  the  basis  for 
the  further  evaluation  of  the  chemistry  of  n-propyl  benzene  and  1 -methyl  naphthalene  as 
described  in  subsequent  sections. 

There  is  a  comparatively  rich  literature  on  the  oxidation  of  toluene  as  compared  to  other 
aromatic  molecules.  Emdee  et  al.  [39]  developed  a  model  for  high  temperature  toluene 
oxidation  utilizing  previous  benzene  and  toluene  schemes  and  validated  it  with  flow  reactor 
experiments  at  temperatures  from  1100  to  1200  K  at  atmospheric  pressure.  Klotz  et  al.  [65] 
updated  the  Emdee  et  al.  [39]  model  and  validated  it  against  toluene-butane  fuel  blends. 
Dagaut  et  al.  [66]  developed  a  model  that  describes  toluene  oxidation  and  validated  it 
against  jet-stirred  reactor  data  for  a  temperature  range  of  1000  -  1375  K  and  equivalence 
ratios  from  0.5  -  1.5.  Lindstedt  and  Maurice  [67]  developed  a  comprehensive  toluene 
mechanism  and  validated  it  against  experimental  data  from  plug  flow  reactors,  shock  tubes, 
counterflow  diffusion  flames  and  premixed  flames.  Sivaramakrishnan  et  al.  [68]  developed 
a  detailed  chemical  model  based  on  an  earlier  literature  model  for  toluene  oxidation  and 
validated  it  against  shock  tube  experimental  data  at  temperatures  in  the  region  1200  -  1500 
K,  over  a  wide  pressure  range  (25  -  610  bar)  and  for  stoichiometries  of  O  =  1  and  0  =  5. 


18 


Final  Report  on  Grant  FA8655-06-1-3052 


A  number  of  experimental  studies  reported  in  the  literature  describe  toluene  oxidation  or 
thermal  decomposition.  Ignition  delay  times  in  argon  mixtures  were  measured  by  Burcat  et 
al.  [69]  for  toluene  concentration  of  0.5  to  1.5%,  4.48  to  13.45%  oxygen,  temperatures  in 
the  range  of  1339  -  1797  K  and  with  reflected  shock  pressures  from  1.95  to  8.85  atm.  More 
recent  studies  on  toluene  oxidation  were  performed  by  Vasudevan  et  al.  [70]  and  Davidson 
and  co-workers  [71].  Vasudevan  et  al.  [70]  measured  ignition  delay  times  and  OH  radical 
concentrations  in  toluene/02/Ar  mixtures  behind  reflected  shock  waves  at  high  temperatures 
(1400  -  2000  K)  and  low  pressures  (1-4  atm)  for  equivalence  ratios  of  0.5  -  1.875  with 
toluene  concentrations  of  0.025  -  0.5%.  The  measured  toluene  ignition  delay  times  were 
compared  with  the  models  of  Pitz  et  al.  [72],  Dagaut  et  al.  [66]  and  Lindstedt  and  Maurice 
[67].  At  low  to  moderate  temperatures,  the  data  agrees  with  the  model  of  Dagaut  et  al.  [66] 
and  Pitz  et  al.  [72]  whereas  at  higher  temperatures  the  model  by  Dagaut  et  al.  [66]  shows  a 
better  agreement  with  measurements.  The  measurements  of  ignition  delay  times  against  fuel 
concentration  shows  that  the  ignition  time  falls  as  the  fuel  concentration  increases  and  this 
trend  is  also  supported  by  Burcat  et  al.  [69].  The  Lindstedt  and  Maurice  [67]  model 
captures  the  trend  of  the  experimental  model  but  with  lower  ignition  delay  times.  The  Pitz 
et  al.  [72]  and  Dagaut  et  al.  [66]  models  can  predict  the  ignition  delays  at  low  fuel 
concentrations  but  show  a  significant  disagreement  at  high  fuel  concentrations. 

Davidson  et  al.  [71]  measured  OH  concentrations  and  ignition  delay  times  in  a  shock  tube 
for  toluene/air  at  low  temperatures  (855  -  1269  K)  and  high  pressures  (14  -  59  atm)  for 
equivalence  ratios  of  0.5  and  1  in  synthetic  air.  The  ignition  delays  were  compared  against 
the  models  of  Pitz  et  al.  [72]  and  Dagaut  et  al.  [66].  The  model  of  Dagaut  et  al.  [66] 
predicts  ten  times  longer  ignition  delays  than  the  experiments  of  Davidson  et  al.  [71]  and 
the  model  of  Pitz  et  al.  [72]  predicts  ignition  delay  times  up  to  two  times  longer  than  the 
measurements.  However,  for  O  =  0.5  the  model  of  Pitz  et  al.  [72]  captures  the  data  well. 

The  toluene  initiation  reactions  have  been  the  target  of  several  studies.  Pamidimukkala  et 
al.  [73]  performed  shock  tube  studies  with  time-of-flight  mass  spectrometry  and  laser 
schlieren  densitometry  for  a  temperature  range  of  1550  -  2000  K  and  pressures  of  0.2  -  0.5 
atm.  Their  study  supports  the  dominance  of  C-C  scission  in  toluene  dissociation  leading  to 
the  production  of  methyl  and  phenyl  radicals  (39)  against  the  C-H  scission  leading  to  the 
benzyl  radical  and  the  hydrogen  atom  (38). 


C7H8 

&1079 

*  c7H7 

+  H 

(38) 

c7h8 

^1 080 

-  c6H5 

+  ch3 

(39) 
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In  addition,  it  was  suggested  that  the  stability  of  the  benzyl  radical  leads  to  a  partial 
equilibrium,  whereas  the  less  stable  phenyl  radical  makes  reaction  (39)  irreversible  and 
hence  a  major  route  in  the  thermal  dissociation  of  toluene.  Colket  and  Seery  [62]  performed 
shock  tube  studies  fortemperature  range  of  1100  -  2700  K  and  suggested  that  (39)  always 
controls  the  distribution  of  products  from  toluene  dissociation. 

Rao  and  Skinner  [74]  monitored  the  formation  of  hydrogen  atoms  under  pyrolytic 
conditions  behind  shock  waves  in  argon  mixtures  at  a  total  pressure  of  0.4  atm  and  for  a 
temperature  range  of  1450  -  1790  K.  Rate  constants  for  the  two  major  toluene  dissociation 
channels  were  proposed  and  it  was  concluded  that  the  hydrogen  abstraction  reaction  (38)  is 
the  only  important  initiation  step.  The  study  was  also  supported  by  Muller-Markgraf  and 
Troe  [75],  who  studied  the  benzyl  radical  absorption  during  its  decay  behind  shock  waves. 
Added  to  this,  Brouwer  et  al.  [76],  utilising  UV  spectroscopy  to  study  the  toluene  thermal 
decomposition,  highlighted  the  dominance  of  the  hydrogen  dissociation  channel. 

Braun-Unkhoff  [77]  monitored  the  formation  of  the  hydrogen  atom  during  thermal 
decomposition  in  the  temperature  range  of  1300-1800  K,  pressures  of  1.5  -  7.8  bar  in  argon 
mixtures  and  highlighted  the  importance  of  CH3  formation  (39)  for  the  product  distribution. 
The  results  also  suggested  that  the  hydrogen  abstraction  from  toluene  leading  to  the 
formation  of  benzyl  radical  is  important. 

Rate  constants  were  proposed  in  the  above  studies  for  the  two  dissociation  channels.  For  the 
C-H  bond  scission,  Rao  and  Skinner  [74]  assigned  an  energy  barrier  of  360  kJ/mole,  Braun- 
Unkhoff  et  al  [77]  assigned  a  value  of  374  kJ/mole  and  Brouwer  et  al.  [76]  a  value  of  369 
kJ/mole.  Later  shock  tube  studies  by  Hippier  and  Troe  [78]  revised  the  C-H  bond  scission 
dissociation  channel  and  assigned  a  new  rate  with  a  barrier  of  360  kJ/mole  over  a 
temperature  range  of  1200-1500  K  and  a  value  of  356  kJ/mole  in  a  flow  system  over  a 
temperature  range  of  913-1143  K. 

Eng  et  al.  [79]  investigated  the  thermal  decomposition  of  toluene  at  temperatures  in  the 
range  1350  -  1900  K  and  for  pressures  from  0.1  to  2.0  bar.  It  was  suggested  that  toluene 
dissociation  proceeds  mainly  via  (38)  for  which  a  rate  constant  was  determined  based  on  H- 
atom  detection  via  calibrated  atomic  resonance  absorption  spectroscopy.  The  reaction 
showed  no  pressure  dependence  at  lower  temperatures,  whereas  at  the  highest  temperature  a 
slight  pressure  dependence  was  observed.  Moreover,  Eng  et  al.  [79]  suggested  that  the 
branching  ratio  k38  /  (k38  +  k39)  for  toluene  decomposition  is  affected  by  temperature  and 
pressure. 

A  more  recent  study  of  the  thermal  decomposition  of  toluene  was  performed  by 
Oehlschlaeger  et  al.  [80],  who  investigated  the  contribution  of  the  two  dissociation  channels 
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in  shock  wave  experiments  over  a  temperature  range  of  1400  -  1780  K  and  at  a  pressure  of 
1.5  bar.  The  benzyl  radical  absorption  at  266  nm  was  monitored  during  toluene 
decomposition  in  argon  mixtures  and  the  rates  for  the  two  channels  were  determined. 
Moreover,  it  was  suggested  that  the  above  branching  ratio  between  the  two  channels  varied 
from  0.8  at  1450  K  to  0.6  at  1800  K.  The  findings  were  compared  to  previous  studies  and  it 
was  found  that  the  overall  decomposition  rate  agrees  with  a  deviation  of  less  than  30%  from 
the  results  of  Braun-Unkhoff  et  al.  [77]  and  Eng  et  al.  [79]  and,  hence,  is  faster  than 
proposed  by  Pamidimukkala  et  al.  [73]  and  Rao  and  Skinner  [74].  The  branching  ratio  is  in 
good  agreement  with  Eng  et  al.  [79],  but  in  disagreement  with  Pamidimukkala  et  al.  [73] 
and  Braun-Unkhoff  et  al.  [77]. 

In  the  present  study,  the  comprehensively  validated  detailed  mechanism  developed  by 
Lindstedt  and  Maurice  [67],  with  subsequent  developments  by  Potter  [60],  was  used  as  a 
starting  point  for  the  re-evaluation  of  the  toluene  chemistry  in  light  of  recent  studies.  The 
toluene  sub-mechanism  features  103  elementary  reactions  and  23  species.  Further  analysis 
was  performed  and  improvements  made  by  considering  new  reaction  channels  and  rates. 

3.1  The  Current  Contribution 

The  complete  mechanism  was  validated  and  analysed  against  data  and  conditions 
corresponding  to  i)  shock  tube  data  for  toluene  pyrolysis  [77]  and  ii)  shock  tube  data  for 
toluene  oxidation  [69,  70].  The  rates  of  toluene  pyrolysis  were  analysed  at  conditions 
corresponding  to  the  shock  tube  experiments  of  Braun-Unkhoff  et  al.  [77].  Moreover,  the 
toluene  sub-mechanism  was  analysed  at  conditions  corresponding  to  toluene  oxidation  in 
argon  mixtures  behind  reflected  shock  waves  (e.g.  Vasudevan  et  al.  [70]  and  Burcat  et  al. 
[69]).  The  validation  and  analysis  of  the  current  model  was  performed  in  order  to  analyse 
the  principal  reaction  paths  and  to  evaluate  the  new  reaction  steps  and  reaction  rate  updates. 
All  the  elementary  reaction  steps  are  assumed  reversible  with  the  reverse  rates  computed 
via  chemical  equilibrium.  The  following  decomposition  pathways  were  updated  using  the 
reaction  rates  determined  by  Oehlschlaeger  et  al.  [80]. 
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=  C7H7  +  H 

(38) 

C7H8 

=  C6H5  +  ch3 

(39) 

CvH7 

=  C7H6  +  H 

(37) 

c7h8 

+ 

H 

=  C7H7  +  h2 

(40) 

c7h8 

+ 

H 

=  C6H6  +  ch3 

(41) 

c7h8 

+ 

ch3 

=  c7h7  +  ch4 

(44) 

Final  Report  on  Grant  FA8655-06-1-3052 


21 


Another  addition  involves  the  thermal  decomposition  of  the  benzylperoxy  radical  (50).  The 
reaction  rate  was  adjusted  to  the  reverse  rate  of  the  corresponding  dissociation  reaction  of 
the  phenylperoxy  radical  (31)  from  DiNaro  et  al.  [81]. 

C7H700  =  C7H70  +  O  (50) 

C6H500  =  C6H50  +  O  (-31) 

Reaction  (50)  was  tested  under  the  shock  tube  conditions  of  Davidson  et  al.  [71]  and  was 
found  to  play  a  significant  role  in  the  evaluation  of  the  OH  profile  and  reduced  the  ignition 
delay  time  by  approximately  22%  as  compared  to  the  starting  mechanism.  The  current 
model  still  fails  to  accurately  capture  low  temperature  (<  1269  K)  ignition  data  and  the 
following  reaction  path  was  also  considered. 


C7H7  +  02 

=  C7H700 

(36) 

C7H700 

=  C7H6OOH 

(Rl) 

C7H6OOH 

=  C7H60  +  OH 

(R2) 

Reaction  (Rl)  was  found  to  be  the  rate-limiting  step  of  the  reaction  sequence.  Due  to  lack 
of  accurate  pathway  information,  the  global  step  given  by  reaction  (49)  was  added  to  the 
mechanism  with  a  rate  and  barrier  based  on  suggestions  of  Hunter  et  al  [82],  with  the 
relative  isomerization  C2H500  =  C2H4OOH  used  as  a  reference  for  estimation  purposes. 

C7H700  =  C7H60  +  OH  (49) 

The  three-step  reaction  pathway  (36,  Rl,  R2)  follows  the  study  of  Zellner  and  Ewig  [83], 
Walch  [84]  and  Clothier  et  al.  [85]  for  the  arguably  related  sequence  (R3  to  R6) 


CH3  +  02 

=  CH3OO 

(R3) 

CH3OO 

=  CH30  +  0 

(R4) 

CH3OO 

=  CH2OOH 

(R5) 

CH2OOH 

=  CH2O  +  OH 

(R6) 

The  reaction  of  methyl  with  molecular  oxygen  producing  the  methylperoxy  radical  leads 
either,  via  decomposition,  to  the  methoxy  radical  and  atomic  oxygen  (R4)  or,  via 
isomerization,  to  CH2OOH  (R5).  The  later  step  is  potentially  important  pathway  under 
radical  depleted  conditions  (e.g.  during  ignition).  The  activation  energy  for  the  transition 
state  (R4)  was  found  to  be  243  kJ/mol  whereas  the  activation  energy  for  the  transition  state 
(R5)  was  only  160  kJ/mol.  Hence,  it  was  suggested  that  the  production  of  CH20  +  OH  via 
(R6)  occurs  with  a  smaller  activation  barrier  and  dominates  at  temperatures  below  2800  K. 
The  two  different  channels  for  the  CH300  decomposition  come  from  different  electronic 
states.  The  species  can  exist  either  in  the  X3Xg-  ground  state  (R4),  where  the  closest  2p 
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orbital  of  the  oxygen  is  double  occupied  from  the  electrons  and  1,3-hydrogen  migration  is 
not  possible,  or  in  the  lAg  state  (R5)  where  the  same  2p  orbital  of  the  oxygen  is  singly 
occupied  and  hydrogen  migration  is  favourable  [84].  A  schematic  representation  of  the  the 
ground  and  excited  state  of  the  oxygen  atom  on  the  benzylperoxy  radical  is  shown  in  Fig. 
3.1.  However,  according  to  Clothier  et  al.  [85]  the  activation  energy  of  reactions  (R3-R4)  is 
approximately  121  kJ/mol  less  that  the  value  estimated  by  Zellner  and  Ewig  [83]. 


Figure  3.1  Schematic  representation  of  electronic  states  of  the  two  benzoperoxy  radicals. 

It  is  obvious  that  the  above  analogy  is  in  need  of  refinement  and  that  the  current  rates  are 
subject  to  uncertainties.  However,  the  pathway  increases  in  significance  under  lower 
temperature  conditions  and  the  current  study  shows  that  further  work  is  desirable.  Another 
reaction  step  (2)  that  plays  a  significant  role  in  both  toluene  pyrolysis  and  oxidation 
processes  was  added  with  a  rate  adopted  from  Colket  and  Seery  [62]. 

C5H5  +  C2H2  =  C7H7  (2) 

A  rate  from  Muller-Markgraf  and  Troe  [75]  was  also  tested  for  the  reverse  reaction  and  was 
found  to  lead  to  an  increase  in  ignition  delays  of  approximately  20%. 

The  thermal  decomposition  of  toluene  in  shock  tubes  was  studied  at  conditions 
corresponding  to  experiments  performed  by  Braun-Unkhoff  et  al.  [77]  in  order  to  evaluate 
the  branching  ratio  of  the  main  toluene  decomposition  reaction  channels  and  to  monitor  the 
temporal  evolution  of  the  hydrogen  radical.  Braun-Unkhoff  et  al.  [77]  studied  the  thermal 
decomposition  in  argon  mixtures  and  measured  H  profiles  at  fuel  concentrations  of  2  -  19.3 
ppm,  temperatures  of  1515  -  1655  K  and  pressures  of  1.89  -  1.93  bar  (Table  3.1).  An 
attempt  was  made  to  model  the  experimental  profiles  using  a  ten-step  reaction  scheme  to 
which  adjustments  had  to  be  made  in  order  to  achieve  the  appropriate  fit  to  the  profiles. 

The  analysis  outlined  below  follows  Lindstedt  and  Maurice  [67]  and  was  performed  to 
assess  the  consistency  of  the  revised  principal  reaction  paths  mentioned  above. 
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C7H8  (ppm) 

T  [K] 

P  [bar] 

1 

2.0 

1515 

1.85 

2 

2.8 

1655 

1.89 

3 

3.0 

1585 

1.93 

4 

19.3 

1555 

1.92 

Table  3.1  Experimental  and  modelling  conditions  for  toluene  thermal  decomposition 

The  computed  hydrogen  profiles  for  all  the  experimental  conditions  are  shown  in  Figs.  3.2 
-3.5.  The  modelled  hydrogen  profiles  are  reasonably  well  captured  and  show  that  the 
adopted  branching  ratio  and  rates  determined  by  Oehlschlaeger  et  al.  [80]  are  satisfactory. 


Figure  3.2  Hydrogen  radical  concentration  against  time  for  toluene  thermal  decomposition  with  an 
initial  fuel  concentration  of  2.0  ppm,  T  =  1515  K  and  P  =  1.85  bar.  The  solid  line  indicates  the 
current  computations  and  the  circles  indicate  the  experimental  data  from  Braun-Unkhoff  et  al.  [77]. 

The  data  used  to  validate  the  toluene  mechanism  under  oxidation  conditions  was  obtained 
from  Vasudevan  et  al.  [70],  who  monitored  the  time  resolved  OH  radical  concentration 
profiles  in  toluene/02/Ar  mixtures  in  shock  tubes.  Vasudevan  et  al.  [70]  measured  OH 
concentration  profiles  over  a  wide  range  of  conditions.  The  study  provides  unique  time- 
resolved  data  of  a  complementary  nature  to  the  pyrolysis  experiments  discussed  above.  In 
the  present  work,  computations  were  performed  corresponding  to  data  with  a  temperature  of 
1689  K  and  pressure  of  1.79  atm  (®  =  1,  0.1%  C7H8,  0.9%  02,  99%  Ar).  The  computed 
profile  is  shown  in  Fig.  3.6  along  with  the  experimental  data  of  Vasudevan  et  al.  [70]. 

There  is  excellent  agreement  between  the  measurements  and  the  current  toluene  model.  The 
OH  profile  can  be  divided  in  three  regions;  The  first  shows  the  rapid  increase  in  the  OH 
concentration  due  to  toluene  decomposition,  the  second  where  OH  rises,  indicating  the 
presence  and  impact  of  the  chain  branching  reactions  and  the  third  where  the  OH 
production  is  close  to  zero  [70]. 
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Figure  3.3  Hydrogen  radical  concentration  against  time  for  toluene  thermal  decomposition  with  an 
initial  fuel  concentration  of  2.8  ppm,  T  =  1655  K  and  P  =  1.89  bar.  The  solid  line  indicates  the 
current  computations  and  the  circles  indicate  the  experimental  data  from  Braun-Unkhoff  et  al.  [77]. 


Figure  3.4  Hydrogen  radical  concentration  against  time  for  toluene  thermal  decomposition  with  an 
initial  fuel  concentration  of  3.0  ppm,  T  =  1585  K  and  P  =  1.93  bar.  The  solid  line  indicates  the 
current  computations  and  the  circles  indicate  the  experimental  data  from  Braun-Unkhoff  et  al.  [77]. 


Figure  3.5  Hydrogen  radical  concentration  against  time  for  toluene  thermal  decomposition  with  an 
initial  fuel  concentration  of  19.3  ppm,  T  =  1555  K  and  P  =  1.92  bar.  The  solid  line  indicates  the 
current  computations  and  the  circles  indicate  the  experimental  data  from  Braun-Unkhoff  et  al.  [77]. 
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A  sensitivity  analysis  was  also  performed  and  the  current  model  shows  that  the  OH  profile 
is  sensitive  to  the  chain  branching  reaction  O  +  H2  — >  OH  +  O.  The  rate  of  Sun  et  al.  [86], 
adopted  from  Baulch  et  al.  [87],  was  tested  for  this  reaction  step  and  found  to  delay  the 
onset  of  ignition  by  approximately  17%  compared  to  the  rate  of  Li  et  al  [88]  adopted  from 
Sutherland  et  al  [89].  A  further  analysis  has  been  presented  by  Gkagkas  and  Lindstedt  [90]. 
It  was  also  shown  that  the  reaction  step  is  responsible  for  12%  of  the  OH  radical 
production.  In  addition,  the  chain  branching  step  H  +  02  — >  OH  +  O  is  responsible  for  67% 
of  the  OH  production  and  the  rate  has  been  studied  extensively  in  the  past  (e.g.  [91-93]). 
Reaction  (2)  also  plays  significant  role  in  the  OH  evolution  as  shown  in  Fig.  3.7.  The 
absence  of  this  reaction  step  causes  a  delay  in  the  OH  ignition  by  around  20%. 

C5H5  +  C2H2  =  C7H7  (2) 

A  further  example  of  a  OH  concentration  profile  versus  time  is  shown  in  Fig.  3.8  for  a 
stoichimetric  250  ppm  toluene  mixture  at  T  =  1648  K  and  P  =  2.03  atm. 

Ignition  delay  times  were  measured  by  Vasudevan  et  al  [70]  in  toluene/02/Ar  mixtures  for 
a  temperature  range  1510  -  1818  K  with  0  =  1  (0.1%  toluene,  0.9%  oxygen,  99%  argon) 
and  a  pressure  of  1  atm  and  defined  the  ignition  delay  time  as  the  time  needed  for  the  OH 
radical  concentration  to  reach  50%  of  the  peak  value  by  setting  as  zero  the  time  of  the 
arrival  of  the  reflected  shock.  Ignition  delay  times  were  also  measured  for  different  fuel 
mole  fractions  varying  from  2xl0"4  -  2xl0"2  with  O  =  1,  at  a  temperature  of  1600  K  and  a 
pressure  of  1  atm.  Variations  of  the  ignition  delay  times  with  temperature  and  fuel 
concentrations  are  presented  in  Figs.  3.9  -  3.10.  As  is  evident  in  Fig.  3.9,  the  current 
toluene  model  agrees  well  with  measurements  at  higher  temperatures. 


Figure  3.6  Time  resolved  OH  radical  concentrations  obtained  in  Toluene/Oxygen/ Argon  mixtures  in 
a  shock  tube  with  ®  =  1  (0.1%  toluene,  0.9%  oxygen),  T  =  1689  K  and  P  =  1.79  atm.  The  circles 
indicate  the  measurements  from  Vasudevan  et  al  [70]  and  the  solid  line  indicates  the  scurrent 
simulation. 
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Figure  3.7  Sensitivity  analysis  of  time  resolved  OH  radical  concentrations  obtained  in 
Toluene/Oxygen/ Argon  mixtures  in  shock  tube  with  O  =  1  (0.1%  toluene,  0.9%  oxygen),  T  =  1689  K 
and  P  =  1.79  atm.  The  circles  indicate  the  measurements  from  Vasudevan  et  al.  [70],  the  solid  line 
indicates  the  current  model,  the  dashed  line  indicates  the  current  model  with  the  rate  of  Sun  et  al 
[86]  for  the  reaction  O  +  H2  — >  OH  +  O  and  the  dashed  dotted  line  indicates  the  current  model  in  the 
absence  of  the  C5H5  +  C2H2  — >  C7H7  reaction. 


The  dependence  of  ignition  delay  times  on  fuel  concentration  is  shown  in  Fig.  3.10. 
Ignition  delay  times  were  also  measured  by  Burcat  et  al  [69]  for  toluene/oxy gen/ argon 
mixtures  for  a  temperature  range  of  1379  -  1785  K,  pressures  of  1.96  -  2.81  atm  for  a 
stoichiometric  mixture  (0.497%  C7H8,  4.48%  02,  95.023%  Ar).  The  comparison  between 
the  model  and  the  experimental  data  is  shown  in  Fig.  3.11. 


Figure  3.8  Time  resolved  OH  radical  concentrations  obtained  in  Toluene/Oxy gen/ Argon  mixtures  in 
a  shock  tube  with  0  =  1  (0.025%  toluene,  0.225%  oxygen),  T  =  1648  K  and  P  =  2.03  atm.  The 
circles  indicate  the  measurements  from  Vasudevan  et  al  [70]  and  the  solid  line  indicates  the  current 
simulation. 
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Figure  3.9  Ignition  delay  times  of  toluene  with  O  =  1  (0.1%  toluene,  0.9%  oxygen,  99%  argon)  and  P 
=  1  atm.  Symbols  indicate  measurements  from  Vasudevan  et  al.  [70]  and  the  solid  line  indicates  the 
current  computations. 


Figure  3.10  Ignition  delay  times  for  toluene  mole  fractions  of  2x1  O'4  -  lxl O'2  with  0=1 
(toluene/oxygen/argon),  T  =  1600  K  and  P  =  1  atm.  The  solid  squares  indicate  the  measurements 
from  Vasudevan  et  al.  [70],  the  open  circles  indicate  experimental  data  from  Burcat  et  al  [69]  and 
the  solid  line  indicates  the  current  computations. 


Figure  3.11  Ignition  delay  times  for  toluene  obtained  in  a  shock  tube  with  0  =  1  (0.497%  C7H8, 
4.48%  02,  95.023%  Ar)  and  P  =  2.28  atm.  The  circles  indicate  experimental  data  from  Burcat  et  al 
[69]  and  the  solid  line  indicates  the  current  computations. 
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3.2  Path  Analysis 

A  reaction  rate  analysis  was  performed  for  the  toluene  decomposition  in  shock  tubes  at  the 
experimental  conditions  of  Braun-Unkhoff  et  al.  [77],  as  discussed  above,  and  the  major 
reaction  pathways  were  determined.  For  the  case  of  a  toluene  concentration  of  2  ppm  at 
1515  K  and  1.85  bar,  the  toluene  decomposition  is  controlled  by  two  major  reaction 
channels.  The  C-H  fission  to  the  benzyl  radical  and  atomic  hydrogen  (38)  is  responsible  for 
67%  of  the  fuel  decomposition.  The  reaction  channel  that  proceeds  via  C-C  scission  to  the 
formation  of  phenyl  and  methyl  radical  accounts  for  27%  of  the  toluene  decomposition.  The 
rates  determined  by  Oehlschlaeger  et  al.  [80]  were  applied  to  these  channels. 

The  current  findings  highlight  the  dominance  of  reaction  (38)  in  agreement  with  studies  of 
Muller-Markgraf  and  Troe  [75],  Brouwer  et  al.  [76],  Rao  and  Skinner  [74]  and  Hippier 
[78].  There  is,  however,  a  disagreement  with  Pamidimukkala  et  al.  [73]  who  proposed 
reaction  (39)  as  the  main  decomposition  channel.  As  shown  in  the  work  by  Oehlschlaeger  et 
al.  [80],  the  overall  decomposition  rate  agrees  with  Braun-Unkhoff  [77],  Eng  et  al.  [79]  and 
Luther  et  al.  [94]  with  a  deviation  of  less  than  30%.  Recommendations  of  Pamidimukkala 
et  al.  [73]  and  Rao  and  Skinner  [74]  are  significantly  slower.  The  benzyl  radical 
decomposes  (90%)  via  reaction  (2)  leading  to  the  cyclopentadienyl  radical  and  acetylene. 

C5H5  +  C2H2  =  C7H7  (2) 

The  cyclopentadienyl  radical  leads  to  the  formation  of  the  propargyl  radical  and  acetylene 
(60%)  or  isomerizes  to  the  linear  form  of  C5H5(L)  (38%) . 

C5H5  =  C3H3  +  C2H2  (3) 

C5H5  =  C5H5(L)  (8) 

Both  decomposition  channels  (3)  and  (8)  were  assigned  rate  constants  computed  from 
potential  energy  surfaces  using  variable  transition  state  theory  and  Rice-Ramsperger- 
Kassel-Marcus/master  equation  approaches  [15].  A  faster  rate  of  Kern  et  al.  [36]  was  also 
tested  but  lead  to  overproduction  of  the  hydrogen  atom  by  20%.  The  other  main  product  of 
the  toluene  thermal  decomposition,  the  phenyl  radical,  undergoes  isomerization  to  the  linear 
C6H5(B)  or  undergoes  hydrogen  fission  reactions  that  follow  a  route  from  C6H5  — >  C6H4  — ► 
C6H4L  — ►  C6H3  — ►  C6H2  and  both  channels  lead  to  the  formation  of  acetylene  via  1,3- 
butadiyne.  A  schematic  representation  of  the  toluene  thermal  decomposition  for  2.0  ppm  of 
fuel  at  1515  K  and  1.85  bar  is  shown  in  Fig.  3.12.  The  temporal  evolution  of  the  hydrogen 
radical  is  controlled  predominantly  by  reaction  (38),  which  accounts  for  92%  of  the 
production. 

C7H8 


C7H7  +H 


(38) 
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The  concentration  of  the  hydrogen  radical  is  also  sensitive  to  reactions  (40),  (17),  (18)  and 
(41)  which  are  responsible  for  27%,  23%,  12%  and  11%  of  the  consumption  respectively. 
The  rates  for  reactions  (40)  and  (41)  were  adopted  from  Oehlschlaeger  et  al.  [80].  The  rate 
for  the  hydrogen  assisted  hydrogen  abstraction  reaction  (17)  was  adopted  from  Leung  and 
Lindstedt  [58].  The  benzene  formation  reaction  via  hydrogen  addition  to  the  phenyl  radical 
(18)  was  assigned  a  rate  proposed  by  Baulch  et  al.  [91]. 


c7h8 

+  H 

=  c7H7 

+  h2 

(40) 

c7h8 

+  H 

=  c6H6 

+  ch3 

(41) 

c6h5 

+  H 

=  c6H4 

+  h2 

(17) 

c6h5 

+  H 

=  c6H6 

(18) 

The  impact  of  reaction  (40)  on  the  hydrogen  profile  was  also  noted  by  Braun-Unkhoff  et  al 
[77],  who  considered  it  in  the  context  of  a  bringing  agreement  between  their  measured  and 
computed  hydrogen  profiles.  The  rates  proposed  by  Braun-Unkhoff  et  al.  [77]  and  Baulch 
et  al.  [91]  were  tested  for  reaction  (40)  and  caused  an  increase  of  the  hydrogen 
concentration  -20%  compared  to  the  rate  proposed  by  Oehlschlaeger  et  al.  [80] 

The  reaction  paths  for  toluene  oxidation  under  shock  tube  conditions  were  analysed  at  the 
conditions  of  Vasudevan  et  al.  [70],  who  measured  the  temporal  evolution  of  the  OH 
radical  for  a  stoichiometric  mixture  of  0.1%  toluene  and  0.9%  oxygen  in  argon  mixture  at 
1689  K  and  1.79  atm.  The  overall  toluene  decomposition  is  mainly  controlled  by  four 
reactions  as  discussed  above.  These  include  the  two  main  thermal  decomposition  pathways 
leading  to  the  formation  of  benzyl  (38)  and  the  phenyl  radical  (39),  which  are  jointly 
responsible  for  around  35%  of  the  fuel  consumption,  and  the  two  hydrogen  assisted 
hydrogen  abstraction  reactions  (40)  and  (41)  leading  to  the  formation  of  benzyl  (29%)  and 
phenyl  radicals  (10%). 

According  to  the  temporal  evolution  of  the  OH  radical  concentration,  the  fuel 
decomposition  follows  three  stages.  The  first  region  where  initiation  reactions  occur  and  the 
fuel  is  decomposed,  the  second  region  where  ignition  occurs  due  to  chain  branching 
reactions  and  the  third  region  with  zero  net  OH  production.  For  the  specific  test  case 
considered,  the  regions  are  separated  for  times  up  to  100  ps,  250  ps  and  800  ps.  Hence,  the 
contribution  of  each  decomposition  step  varies  with  time.  Although  reaction  (38)  plays  a 
significant  role  in  the  overall  fuel  consumption,  it  is  predicted  that  it  is  responsible  for  20% 
in  the  first  stage,  it  has  no  impact  on  the  subsequent  stages  of  the  fuel  oxidation  process. 
Reaction  (39)  is  responsible  for  19%,  12%  and  14%  of  the  fuel  decomposition  over  the 
three  stages. 
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The  contribution  of  reaction  (40)  to  the  fuel  consumption  increases  over  time  from  27%  in 
the  first  region  to  45%  in  the  second  and  third  regions  respectively.  The  reaction  consumes 
the  H  radical  and  leads  to  the  formation  of  the  less  reactive  benzyl  radical  and  molecular 
hydrogen.  It  must  be  noted  that  reactions  that  consume  reactive  radicals  during  fuel 
oxidation  can  inhibit  ignition.  Reaction  (41)  shows  an  increasing  impact  on  the  fuel 
consumption  from  9%  at  the  first  stage  to  15%  up  to  the  end  of  the  fuel  decomposition. 

The  benzyl  radical  produced  from  reactions  (38)  and  (40),  undergoes  thermal  dissociation 
that  leads  to  the  production  of  the  cyclopentadienyl  radical  and  acetylene  via  reaction  (2) 
which  accounts  for  80%  of  the  benzyl  consumption. 

C5H5  +  C2H2  =  C7H7  (2) 

Benzene  is  mostly  produced  (50%)  via  reaction  (41)  and  consumed  (38%)  via  reaction  (25) 
for  which  a  rate  from  Leung  and  Lindstedt  [58]  was  adopted.  A  rate  of  Leung  and  Lindstedt 
was  also  assigned  to  reaction  (22)  which  is  responsible  for  30%  of  the  benzene 
consumption.  The  production  of  the  phenoxy  radical  and  the  hydrogen  atom  via  oxygen 
attack  on  benzene  is  responsible  for  28%  of  the  consumption  of  the  latter  with  the  rate 
adopted  from  DiNaro  et  al.  [81]. 


cvh8 

+  H 

=  c6h6 

+  ch3 

(41) 

c6h6 

+  OH 

=  c6h5 

+  h2o 

(25) 

c6h6 

+  H 

=  c6h5 

+  h2 

(22) 

c6h6 

+  0 

=  c6h5o 

+  H 

(23) 

The  phenyl  radical  is  initially  produced  (75%)  from  the  thermal  decomposition  of  toluene 
via  reaction  (39).  As  most  of  the  fuel  is  consumed,  secondary  reactions  occur  and  the 
phenyl  radical  pool  is  produced  via  reactions  (24,  22)  that  were  assigned  rates  from  Leung 
and  Lindstedt  [58].  The  contribution  of  reaction  (24)  to  the  phenyl  production  remains  the 
same  up  to  800  ps  (47%),  but  reaction  (22)  is  responsible  for  between  38%  and  48%  of  the 


production. 

C7H8 

=  c6h5 

+  ch3 

(39) 

C6H6  +  OH 

=  c6h5 

+  h2o 

(24) 

C6H6  +  H 

=  c6h5 

+  h2 

(22) 

The  consumption  of  the  phenyl  radical  proceeds  mainly  (54%)  via  reaction  (19)  featuring 
molecular  oxygen  attack  on  the  phenyl  radical  leading  to  the  phenoxy  radical  and  an 
oxygen  atom.  A  rate  from  Frank  et  al.  [95]  was  adopted  for  this  step.  In  addition,  a 
hydrogen  atom  is  abstracted  from  the  ring  (18%)  via  hydrogen  attack  (17)  that  leads  to  the 
production  of  C6H4  and  molecular  hydrogen.  The  rate  assigned  to  reaction  (17)  was 
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obtained  from  Leung  and  Lindstedt  [58].  Approximately  12%  of  the  phenyl  radical 
consumption  occurs  via  reaction  (20),  which  involves  molecular  oxygen  addition  to  the 
phenyl  radical  leading  to  the  formation  of  the  phenyl  peroxy  radical. 


Figure  3.12  Predicted  toluene  thermal  decomposition  pathways  corresponding  to  the  study  of  Braun- 
Unkhoff  et  al.  [77]  with  an  initial  fuel  concentration  of  2.0  ppm,  T  =  1515  K  and  P  =  1.85  bar. 


C6H5 

+  02 

=  c6h5o  +  0 

(19) 

c6h5 

+  H 

=  c6h4  +  h2 

(17) 

c6h5 

+  02 

=  c6h5oo 

(20) 

The  phenoxy  radical  is  mainly  produced  (58%)  from  reaction  (19)  and  leads  (83%)  to 
cyclopentadienyl  and  carbon  monoxide  (28)  and  with  a  fraction  (15%)  responsible  for  the 
formation  of  phenol  (29).  The  rate  for  the  carbon  monoxide  expulsion  was  adopted  from 
Leung  and  Lindstedt  [58]  and  the  rate  for  the  phenol  formation  from  DiNaro  et  al.  [81]. 
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C6H50  =  C5H5  +  CO  (28) 

C6HsO  +  H  =  C6H5OH  (29) 

The  consumption  of  C6H4  follows  the  reaction  route  C6H4  — ►  C6H4L  — >  C6H3  — >  C6H2  — ► 
C4H2  —>  C4H20  — ►  C2H2  leading  to  acetylene.  The  phenylperoxy  radical  decays  via  two 
thermal  pathways  (31,  32)  and  leads  to  the  formation  of  reactive  radicals  that  are 
responsible  for  68%  and  22%  of  the  consumption.  Both  decomposition  steps  were  assigned 
the  rates  of  DiNaro  et  al.  [81].  The  p-benzoquinone  formed  via  (32)  is  decomposed  to 
acetylene,  carbon  monoxide  and  formaldehyde  following  the  global  steps  introduced  by  the 
latter  authors. 

C6H500  =  C6H50  +0  (31) 

C6HsOO  =  C6H402  +H  (32) 

The  cyclopentadienyl  radical,  which  is  produced  via  benzyl  radical  decomposition  or  via 
the  phenyl  radical  degradation  route,  is  decomposed  via  four  main  reaction  pathways. 
Approximately  30%  is  thermally  decomposed  leading  to  the  formation  of  acetylene  and  the 
propargyl  radical  and  29%  reacts  with  the  oxygen  atom  leading  to  carbon  monoxide  and 
C4H5(T).  Moreover,  15%  isomerizes  to  a  linear  structure  and  13%  recombines  with 
hydrogen  forming  cyclopentadiene.  Rates  for  reactions  (3),  (4)  and  (8)  were  calculated  by 
Robinson  [15]  from  potential  energy  surfaces  determined  using  variable  transition  state 
theory  and  Rice-Ramsperger-Kassel-Marcus/master  equation  approaches.  A  rate  from  Kern 
et  al.  [36]  was  applied  to  the  hydrogen  recombination  reaction  (13). 


C5H5 

—  C3H3  +  C2H2 

(3) 

C5H5  +  0 

=  C4H5(T)  +CO 

(4) 

C5H5 

=  C5H5(L) 

(8) 

c5h6 

=  C5H5  +  H 

(13) 

Overall,  the  H  atom  profiles  show  generally  good  agreement  with  measurements  under 
pyrolytic  conditions.  Under  oxidation  conditions  and  at  high  temperatures,  the  current 
toluene  chemistry  captures  the  temporal  evolution  of  the  OH  radical  in  excellent  agreement 
with  the  measurements.  Ignition  delay  times  for  stoichiometric  toluene  mixtures  were  also 
computed  and  showed  good  agreement  with  the  experimental  data  sets  of  Vasudevan  et  al. 
(2005)  and  Burcat  et  al.  (1986). 

The  results  outlined  above  are  encouraging  and  suggest  that  the  current  reaction  mechanism 
can  be  applied  to  other  methyl  substituted  aromatic  fuel  compounds  that  form  part  of  real 
and  surrogate  fuel  formulations. 
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4.0  The  Oxidation  of  n-Propylbenzene 

Previous  studies  of  the  oxidation  of  n-alkyl  benzenes  have  highlighted  that  after  side  chain 
removal,  the  fuels  follow  the  same  oxidation  route  as  that  of  benzene.  However,  benzene 
itself  can  not  represent  more  complex  aromatic  fuels  that  are  present  in  commercial  blends 
and  n-propyl  benzene  is  a  potentially  a  good  candidate  representing  the  mono  alkylated  and 
mono  cyclic  aromatic  component  of  fuels  such  as  gasoline,  diesel  and  kerosene  [96].  The 
current  work  was  focused  on  determining  the  side  chain  removal  process  and  the  steps  that 
characterize  it.  Three  routes  were  found:  (i)  Homolysis  -  direct  cleavage  of  the  side  chain 
followed  by  the  oxidation  of  the  remaining  radical,  (ii)  displacement  of  the  alkyl  side  chain 
by  a  radical  species  and  (iii)  abstraction  of  a  hydrogen  atom  from  the  alkyl  group  [97-100]. 
Due  to  the  analogy  between  alkyl  benzenes  and  alkanes  regarding  atomic  hydrogen 
abstractions,  relevant  reactions  in  the  n-propyl  benzene  chemistry  were  based  on  the 
chemistry  of  propane  as  appropriate.  The  current  work  also  evaluates  the  relevance  of  this 
analogy  in  the  context  of  the  prediction  of  the  combustion  behaviour  of  n-propyl  benzene. 

The  n-propyl  benzene  molecule  has  three  primary,  two  secondary  and  two  benzylic 
hydrogen  atoms.  According  to  the  hydrogen  carbon  bonding  rules,  the  benzylic  bond  is  the 
weakest,  as  compared  to  the  primary  and  secondary  bonds,  thus  it  is  easier  to  abstract  a 
hydrogen  atom  from  the  benzyl  carbon  atom.  The  benzylic  bond  strength  is  around  368 
kJ/mol,  the  secondary  397  kJ/mol  and  the  primary  410  kJ/mol.  The  benzylic  C-H  bond 
strength  is  less  than  at  the  other  two  sites  due  to  the  resonance  of  the  resulting  benzyl 
radical.  The  bond  strengths  thus  suggest  that  the  dominant  abstraction  reaction  features  the 
benzylic  hydrogen  atom  [99,  100].  However,  it  is  difficult  for  any  other  oxidation  of  the 
benzylic  radical  to  occur  as  it  is  stereochemically  hindered  [101].  The  active  radical  species 
are  predominantly  H,  O  and  OH  and  the  reactions  may  be  written  as  below,  where  X  =  {X  I 
H,  O,  OH} 


C6H5-CH2-CH2-CH3  +  X  <-►  C6H5-CH-CH2-CH3  +  XH 
C6H5-CH-CH2-CH3  <-►  C6H5-CH=CH2  +  CH3 


(4.1) 

(4.2) 


The  resulting  benzylic  radical  will  undergo  a  beta-scission  to  form  styrene  and  the  methyl 
radical  (4.1)  -  (4.2).  Hence,  the  early  appearance  of  styrene  shows  that  the  above  step  is  a 
significant  fuel  breakdown  route.  Styrene  can  be  further  decomposed  to  C6H5CHCH, 
benzene  or  react  with  oxygen  leading  to  the  formation  of  oxygenated  species  [99]. 


C6H5-CH=CH2  +X  <-►  C6H5-CH=CH  +  XH 
C6H5-CH=CH  <-►  C6H5-C=CH  +  H 


C6H5-CH=CH2  +  H  <-►  C6H6  +ch=ch2 


(4.3) 

(4.4) 

(4.5) 
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C6H5-CH=CH2  +  O  <-►  C6H5-CH2-CHO  (4.6) 

C6H5-CH2-CHO  <->  C6H5-CH2  +  HCO  (4.7) 

The  primary  hydrogen  abstraction  is  as  follows: 

C6H5-CH2-CH2-CH3  +  X  <-►  C6H5-CH2-CH2-CH2  +  XH  (4.8) 

C6H5-CH2-CH2-CH2  <-►  C6H5-CH2  +  CH2=CH2  (4.9) 


The  resulting  benzyl  radical  may  lead  to  (i)  toluene  formation  via  hydrogen  addition,  (ii) 
benzaldehyde  formation  via  oxygen  addition  or  (iii)  side  chain  alkyl  radical  addition. 
Ethylene  appears  in  large  quantities  in  experimental  studies  of  n-propyl  benzene  oxidation 
showing  that  the  primary  hydrogen  atom  abstraction  is  a  significant  path  in  the  fuel 
oxidation  process  [100].  Significant  quantities  of  products  from  primary  and  benzylic 
hydrogen  abstraction  are  measured  and  it  is  expected  that  secondary  hydrogen  abstractions 
occur  due  to  the  fact  that  the  carbon  bond  strength  is  between  the  two  other  cases. 


C6H5-CH2-CH2-CH3  +  X  <-►  C6H5-CH2-CH-CH3  (4.10) 

C6H5-CH2-CH-CH3  <-►  C6H5-CH=CH-CH3  +  H  (4.11) 

C6H5-CH2-CH-CH3  ^C6H5  +  CH2=CH-CH3  (4.12) 

C6H5-CH2-CH-CH3  <-*  C6H5-C(CH2)-CH3  (4.13) 

C6H5-C(CH2)-CH3  <->  C6H5-CH=CH2  +  CH3  (4.14) 


Apart  from  the  two  products  of  reactions  (4.11)  -  (4.12),  isomerization  (4.13)  or  phenyl 
shift  (4.14)  may  also  occur.  Hydrogen  shift  from  the  benzylic  to  the  secondary  carbon  atom 
is  possible  leading  to  the  a-phenyl  propyl  radical  and,  following  oxidation,  to  styrene. 
Instead  of  a  hydrogen  shift,  it  is  also  possible  for  the  phenyl  group  to  be  transferred  to  the 
beta-carbon  atom  (4.14)  [102].  The  product  of  the  phenyl  shift  (4.14)  will  lead  to  styrene 
production.  The  displacement  of  the  propyl  group  by  a  radical  species  (e.g.  a  hydrogen 
atom)  produces  benzene  and  the  propyl  radical  and  the  detection  of  propylene  in  the 
products  provides  evidence  that  the  displacement  reaction  occurs. 


C6H5-CH2-CH2-CH3  +  H  <-►  C6H6  +  C3H7  (4.15) 

The  homolysis  route  involves  carbon-carbon  cleavage  vi  three  possible  routes. 

C6H5-CH2-CH2-CH3  <->  C6H5  +CH2-CH2-CH3  (4.16) 

C6H5-CH2-CH2-CH3  <-►  C6H5-CH2  +  CH2-CH3  (4.17) 

C6H5-CH2-CH2-CH3  <->  C6H5-CH2-CH2  +CH3  (4.18) 

The  homolysis  route  is  known  to  have  large  activation  energies,  hence  it  will  only  be  an 


important  path  at  high  temperatures.  The  benzylic  C-C  bond  is  weaker  than  the  other  two  in 
the  chain  and  its  breakage  is  expected  to  occur  faster. 
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In  this  study,  the  oxidation  steps  of  n-propyl  benzene  were  analysed  and  the  fuel  breakdown 
was  predicted  for  two  different  pressures  and  different  stoichiometries.  The  reaction 
mechanism  was  derived  from  analogies  with  propane  and  toluene  chemistries  and 
validation  was  performed  with  data  obtained  in  shock  tubes  and  jet  stirred  reactors. 


4.1  The  Current  Contribution 

The  updated  n-propyl  benzene  reaction  mechanism  used  here  consists  of  1683  reversible 
reactions  involving  269  species.  The  thermochemical  data  were  obtained  from  literature 
sources  [56]  and,  when  not  available,  were  calculated  with  quantum  mechanical  methods 
using  Gaussian-03  [15]. 

The  oxidation  of  n-propyl  benzene  was  studied  under  Jet  Stirred  Reactor  (JSR)  conditions 
at  atmospheric  pressure  and  predictions  compared  with  measurements  obtained  by  Dagaut 
et  al.  [96].  Three  different  stoichiometries  (®  =  0.5,  1.0  and  1.5)  were  analyzed  for  a 
temperature  range  of  900  -  1250  K  (Table  4.1).  Concentration  profiles  of  all  major  species 
were  computed  at  a  mean  residence  time  of  x  =  70  ms  and  significant  pathways  for  the  fuel 
breakdown  were  identified.  Comparisons  between  the  simulations  and  experimental  data 
sets  from  Dagaut  et  al.  [96]  are  shown  in  Figs.  4.1  -  4.9. 


0 

P  (atm) 

T  (K) 

X02 

Xc9H12 

0.5 

1.0 

900-1200 

0.024 

0.001 

1.0 

1.0 

950-1250 

0.012 

0.001 

1.5 

1.0 

950-1250 

0.008 

0.001 

Table  4.1  Experimental  and  modelling  conditions  for  the  oxidation  of  ^-propyl  benzene  in  a  jet- 
stirred  reactor  at  P  =  1  atm.  The  species  concentrations  correspond  to  mole  fractions 

The  major  species  measured  by  Dagaut  et  al.  [96],  apart  from  02,  CO  and  C02,  were  ethyl 
benzene  (C8H10),  styrene  (C8H8),  toluene  (C7H8),  benzene  (C6H6),  acetylene  (C2H2), 
ethylene  (C2H4),  methane  (CH4)  and  formaldehyde  (CH20). 

The  fuel  decay  is  reasonably  well  captured  for  all  the  three  tested  equivalence  ratios  and,  as 
the  stoichiometry  increases,  a  better  agreement  between  the  model  and  the  measurements  is 
achieved  for  all  the  intermediate  major  species.  The  agreement  between  predictions  and 
measurements  suggests  that  the  n-propyl  benzene  submechanism  is  adequate. 
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Figure  4.1  Concentration  profiles  of  intermediate  species  during  n-propyl  benzene  oxidation  in  a  jet- 
stirred  reactor  with  ®  =  0.5,  P  =  1  atm,  T  =  900  -  1200  K.  The  circles  are  measurements  [96]  and  the 
solid  line  the  current  simulation. 


Figure  4.2  Concentration  profiles  of  intermediate  species  during  n-propyl  benzene  oxidation  in  a  jet- 
stirred  reactor  with  ®  =  0.5,  P  =1  atm,  T  =  900  -  1200  K.  The  circles  are  measurements  [96]  and  the 
solid  line  the  current  simulation. 


Figure  4.3  Concentration  profiles  of  intermediate  species  during  n-propyl  benzene  oxidation  in  a  jet- 
stirred  reactor  with  ®  =  0.5,  P  =1  atm,  T  =  900  -  1200  K.  The  circles  are  measurements  [96]  and  the 
solid  line  the  current  simulation. 
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Figure  4.4  Concentration  profiles  of  intermediate  species  during  n-propyl  benzene  oxidation  in  a  jet- 
stirred  reactor  with  ®  =  1.0,  P  =  1  atm,  T  =  950  -  1250  K.  The  circles  are  measurements  [96]  and  the 
solid  line  the  current  simulation. 


Figure  4.5  Concentration  profiles  of  intermediate  species  during  n-propyl  benzene  oxidation  in  a  jet- 
stirred  reactor  with  ®  =  1.0,  P  =  1  atm,  T  =  950  -  1250  K.  The  circles  are  measurements  [96]  and  the 
solid  line  the  current  simulation. 


Figure  4.6  Concentration  profiles  of  intermediate  species  during  n-propyl  benzene  oxidation  in  a  jet- 
stirred  reactor  with  O  =  1.0,  P  =  1  atm,  T  =  950  -  1250  K.  The  circles  are  measurements  [96]  and  the 
solid  line  the  current  simulation. 
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Figure  4.7  Concentration  profiles  of  intermediate  species  during  n-propyl  benzene  oxidation  in  a  jet- 
stirred  reactor  with  <X>  =  1.5,  P  =1  atm,  T  =  950  -  1250  K.  The  circles  are  measurements  [96]  and  the 
solid  line  the  current  simulation. 


Figure  4.8  Concentration  profiles  of  intermediate  species  during  n-propyl  benzene  oxidation  in  a  jet- 
stirred  reactor  with  ®  =  1.5,  P  =1  atm,  T  =  950  -  1250  K.  The  circles  are  measurements  [96]  and  the 
solid  line  the  current  simulation 


Figure  4.9  Concentration  profiles  of  intermediate  species  during  n-propyl  benzene  oxidation  in  a  jet- 
stirred  reactor  with  =  1.5,  P  =1  atm,  T  =  950  -  1250  K.  The  circles  are  measurements  [96]  and  the 
solid  line  the  current  simulation. 
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A  reaction  rate  analysis  was  performed  for  a  fuel  rich  mixture  (®  =  1.5)  tested  under 
oxidation  conditions  in  jet-stirred  reactor  at  a  temperature  of  1050  K  and  at  atmospheric 
pressure.  Computations  show  that  the  fuel  decomposes  predominantly  (35%)  by  homolytic 
fission  at  the  secondary  carbon  atom  of  the  branch  leading  to  the  formation  of  benzyl  and 
ethyl  radicals.  An  additional  10%  of  the  fuel  is  consumed  via  hydrogen  abstraction  through 
hydrogen  atom  attack  on  the  ‘benzylic’  (primary)  carbon  atom  leading  to  the  1 -propyl 
benzyl  radical  (123).  The  kinetics  of  reaction  (157)  was  estimated  according  to  the 
approach  of  Dean  [13].  The  same  approach  was  applied  to  the  C-C  homolytic  steps  that 
occur  either  at  the  primary  or  at  the  tertiary  carbon  atom  to  estimate  their  kinetics.  The  rate 
for  reaction  (123)  was  based  on  suggestions  of  Dagaut  et  al.  [96].  Approximately  23%  of 
the  total  fuel  consumption  leads  to  the  formation  of  the  1 -propyl  benzyl  radical  via 
hydrogen  atom  abstraction  with  H,  O  and  OH  radicals.  The  primary  homolytic  reaction, 
which  occurs  via  C-C  scission  leading  to  the  phenyl  and  propyl  radicals,  corresponds  to  8% 
of  the  total  fuel  consumption  with  a  rate  assigned  based  on  the  approach  of  Dean  [13]. 
Approximately  16%  of  the  total  fuel  concentration  is  decomposed  via  hydrogen  abstraction 
reactions  leading  to  the  formation  of  the  2-propyl  benzyl  radical  in  comparison  to  8%  of  the 
fuel  that  leads  to  the  3-propyl  benzyl  radical. 

C9H12  =  C7H7  +  C2H5  (157) 

C9H12  +H  =  'CoHn  +  H2  (123) 

The  1 -propyl  benzyl  radical  decomposes  (100%)  to  styrene  via  methyl  radical  abstraction 
with  a  rate  adopted  by  Dagaut  et  al.  [96].  The  step  also  constitutes  the  major  styrene 
formation  pathway.  Once  styrene  was  formed,  five  major  routes  were  detected  that 
contribute  to  its  decay.  Reaction  (59)  was  assigned  a  rate  proposed  by  Maurice  [103]  and 
contributes  21%  to  the  consumption  of  styrene.  The  C8H7(P)  recycles  back  to  styrene 
(78%).  The  formation  of  the  benzyl  radical  and  CHO  (61)  constitutes  the  second  major 
consumption  pathway  (18%)  with  a  reaction  rate  suggested  by  Potter  [60]. 

C8H8  +  OH  =  C8H7(P)  +  H20  (59) 

C8H8  +0  =  C7H7  +  CHO  (61) 

Reactions  (21),  (16)  and  (62)  make  approximately  equal  contributions  to  the  styrene 
consumption  at  13,  12  and  10%  respectively. 


c6h6 

+  c2h3 

II 

n 

00 

a 

00 

+  H 

(21) 

c6h5 

+  c2h4 

II 

n 

00 

a 

oo 

+  H 

(16) 

c8h8 

+  0 

=  c6H5 

+  CH3CO 

(62) 
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The  major  product  of  the  fuel  decay,  the  benzyl  radical,  recombines  with  a  hydrogen  atom 
leading  to  the  formation  of  toluene.  The  recombination  step  is  the  major  toluene  formation 
channel  (90%)  and  was  assigned  a  rate  from  Oehlschlaeger  et  al.  [80].  Approximately  34% 
of  the  benzyl  radical  reacts  with  the  methyl  radical  and  forms  ethyl  benzene  with  the  rate 
proposed  by  Lindstedt  et  al.  [45].  The  reaction  step  is  the  major  ethyl  benzene  formation 
channel  (97%).  Moreover,  20%  of  the  benzyl  radical  recombines  leading  to  C14H14  with  a 
reaction  rate  obtained  from  Oehlschlaeger  et  al.  [80].  The  major  benzyl  consumption  steps 
are  shown  in  Fig.  4.10  and  the  overall  oxidation  steps  in  Fig.  4.11. 

Ethyl  benzene  is  consumed  via  two  major  steps.  The  displacement  of  the  ethyl  branch  via  a 
hydrogen  atom  is  responsible  for  50%  of  the  ethyl  benzene  consumption.  The  second  major 
channel  occurs  via  hydrogen  atom  abstraction  to  the  formation  of  ethylbenzyl  radical 
(34%),  which  essentially  recycles  back  to  ethyl  benzene  (95%). 

C8H10  +  H  =  C6H6  +  C2H5  (64) 

C8H10  +  H  =  C8H9  +  H2  (63) 

Benzene  is  partly  consumed  via  oxygen  atom  addition  and  hydrogen  abstraction  (23) 
leading  to  C6H50  (58%)  with  a  rate  adopted  from  DiNaro  et  al.  [81]  and  (34%)  to  the 
phenyl  radical  via  OH  attack  (25)  with  a  rate  of  Leung  et  al.  [58].  Oxygen  addition  to  the 
phenyl  radical  leads  to  the  production  of  C6H50  (57%)  and  C6H500  (42%).  The  latter 
species  are  responsible  for  cyclopentadienyl  radical  production  and  contribute  70%  and 
25%  respectively. 

C6H6  +0  =  C6HsO  +H  (23) 

C6H6  +  OH  =  C6H5  +  H20  (25) 

As  mentioned  above,  the  major  fuel  consumption  pathway  produces  benzyl  and  ethyl 
radicals  and  also  constitutes  the  major  ethyl  radical  formation  step  (96%).  The  ethyl  radical 
is  consumed  via  hydrogen  abstraction  (95%)  leading  to  the  formation  of  ethylene  (76%). 
Acetylene  is  also  an  important  specie  produced  during  propyl  benzene  oxidation.  The  main 
channel  is  via  the  thermal  decomposition  of  the  cyclopentadienyl  radical  (41%),  which 
features  an  adjusted  reaction  rate  based  on  that  proposed  by  Kern  et  al.  [36].  An  additional 
28%  of  the  acetylene  production  occurs  via  methyl  abstraction  from  the  C3H5(S)  radical 
formed  through  the  following  reaction  route  C9H12  2C9Hn  C3H6  C3H5(S)  C2H2. 

Methane  is  formed  via  reaction  routes  directly  linked  to  the  fuel  and  involves  hydrogen 
abstraction  via  methyl  radical  attack.  The  hydrogen  abstraction  reactions  lead  to  the 
formation  of  ^Hn  and  3C9Hn  with  each  of  these  contributing  20%  to  the  methane 
production.  It  must  be  noted  that  the  methane  production  channel  that  produces  2C9Hn  is 
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responsible  for  not  more  than  3%  of  the  methane  production.  The  rates  assigned  to  the 
following  reactions  are  based  on  the  kinetics  of  propane  suggested  by  Tsang  [104]. 


c9H12 

+  ch3 

=  'C9H„ 

+  ch4 

(133) 

c9h12 

+  ch3 

=  3C9H„ 

+  ch4 

(136) 

c9h12 

+  ch3 

=  2c9H„ 

+  ch4 

(137) 

The  oxidation  of  n-propyl  benzene  was  also  studied  under  Jet  Stirred  Reactor  (JSR) 
conditions  at  a  pressure  of  10  atm  and  validated  using  the  measurements  obtained  by 
Dagaut  et  al.  [96].  The  stoichiometries  tested  were  0.5,  1.0,  1.5  and  2.0  for  a  temperature 
range  of  900  -  1200  K  at  a  residence  time  of  x  =  0.5  sec.  Comparisons  between  the 
simulations  and  experimental  data  sets  from  Dagaut  et  al.  [96]  are  similar  to  those  discussed 
above  for  atmospheric  pressure  and  are  hence  omitted  from  the  current  report. 


Figure  4.10  Benzyl  radical  consumption  paths  for  O  =  1.5,  T  =  1050  K  and  P  =  1  atm  in  a  jet-stirred 
reactor 

Ignition  delay  times  for  n-propyl  benzene  were  also  computed  and  validated  against 
measurements  obtained  by  Eberius  et  al.  [108]  as  shown  in  Fig.  4.12.  The  data  were 
obtained  in  n-propyl  benzene/02/Ar  mixtures  for  a  temperature  range  of  1400  -  1760  K,  a 
pressure  of  5  atm  and  stoichiometries  from  0.8  -  1.142.  The  ignition  delay  time  was  defined 
by  Eberius  et  al.  [108]  as  the  time  needed  for  the  CH-emission  signal  to  reach  its  peak 
value.  Low  dilution  conditions  were  responsible  for  difficulties  in  evaluating  the  measured 
data  due  to  the  change  in  density  during  the  progress  of  reaction.  Hence,  a  series  of 
experiments  for  fuel/02  mixtures  diluted  in  Argon  were  used.  The  set  of  conditions 
computed  are  shown  in  Table  4.2. 
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Figure  4.11.  n-Propyl  benzene  oxidation  pathways  for  <D  =  1.5,  T  =  1050  K  and  P  =  1  atm  in  a  jet 
stirred  reactor 

A  comparison  of  experimental  data  with  computaitons  is  shown  in  Fig.  4.12  and  arguably 
the  current  n-propyl  benzyl  model  provides  satisfactory  predictions  of  the  ignition  delay 
times  at  intermediate  temperatures. 


[C9H12  ](ppm) 

[02]  (ppm) 

[Ar]  (ppm) 

O 

P  (atm) 

T  (K) 

1 

537 

8057 

991406 

0.800 

5 

1437 

2 

545 

8057 

991398 

0.812 

5 

1416 

3 

537 

8057 

991406 

0.800 

5 

1569 

4 

545 

6289 

993166 

1.040 

5 

1468 

5 

544 

6289 

993167 

1.038 

5 

1631 

6 

541 

6289 

993170 

1.032 

5 

1549 

7 

598 

6284 

993118 

1.142 

5 

1558 

8 

536 

6284 

993180 

1.024 

5 

1555 

9 

542 

6284 

993174 

1.035 

5 

1478 

10 

536 

6284 

993180 

1.024 

5 

1759 

11 

538 

6295 

993167 

1.026 

5 

1448 

12 

537 

6295 

993168 

1.024 

5 

1606 

13 

553 

6295 

993152 

1.054 

5 

1550 

Table  4.2  Experimental  and  modelling  conditions  for  ^-propyl  benzene  ignition  delay  times  in  a 
shock  tube 
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Figure  4.12  Ignition  delay  times  of  ^-propyl  benzene  at  P  =  5  atm,  0.71  <  O  <  1.18.  The  circles  are 
measurements  [108]  and  the  solid  line  the  current  simulation. 

A  reaction  rate  analysis  was  performed  at  1448  K  (Case  11,  Table  4.2)  and  major  reaction 
channels  for  the  fuel  breakdown  were  identified  corresponding  to  Fig.  4.13.  The  fuel  is 
solely  decomposed  via  thermal  breakdown  to  phenyl  (46%)  and  benzyl  (51%)  radicals.  The 
rates  assigned  to  the  thermal  decomposition  steps  (156)  and  (157)  were  based  on  the 
approach  of  Dean  [13]. 

C9H12  =  C6H5  +  C3H7(N)  (156) 

C9H12  =  C7H7  +  C2H5  (157) 

The  benzyl  radical  is  predominantly  (65%)  consumed  via  hydrogen  recombination  to 
toluene.  Under  jet  stirred  reactor  conditions  the  formation  of  toluene  is  responsible  for  35% 
of  the  benzyl  radical  consumption.  Most  of  the  toluene  (47%)  recycles  back  to  the  benzyl 
radical  via  hydrogen  abstraction  reactions  via  H  (40)  and  OH  radical  attack  (43),  while  14% 
produces  benzene  via  C-C  scission  (41)  and  12%  is  oxidized  to  OC7H7  (45).  The  rate 
assigned  to  reaction  (40)  was  adopted  from  Oehlschlaeger  et  al.  [80]  and  the  second  benzyl 
radical  formation  channel  (43)  used  a  rate  proposed  by  Baulch  et  al.  [91].  The  rate  proposed 
by  Oehlschlaeger  et  al.  [80]  was  used  for  reaction  (41)  while  the  toluene  oxidation  channel 
(45)  was  assigned  a  rate  suggested  by  Hoffman  et  al.  [109]. 


c7h8 

+  H 

=  c7h7 

+  h2 

(40) 

c7h8 

+  OH 

=  c7h7 

+  h2o 

(43) 

c7h8 

+  H 

II 

n 

Os 

a 

Os 

+  ch3 

(41) 

c7h8 

+  o 

=  oc7h7 

+  H 

(45) 

Most  of  the  benzene  is  consumed  via  hydrogen  abstraction  with  H  (21%)  and  OH  (42%) 
radicals  leading  to  the  phenyl  radical.  An  additional  32%  of  benzene  is  oxidized  to  phenoxy 
via  reaction  with  atomic  oxygen.  The  phenyl  radical  constitutes  the  second  major  direct 
product  from  the  fuel  consumption.  However,  the  production  comes  predominantly  from 
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the  decomposition  of  benzene.  Once  the  phenyl  radical  is  formed,  71%  is  oxidised  to 
C6H50  via  (19)  and  22%  to  C6H500  via  (20)  respectively. 


C6H5  +  o2 

=  c6h5o  +0 

(19) 

C6H5  +  o2 

=  c6H5oo 

(20) 

The  phenoxy  radical  is  further  decomposed  to  cyclopentadienyl  radical  (80%),  which 
subsequently  is  thermally  decomposed  (73%)  to  the  propargyl  radical  and  acetylene. 
Approximately  34%  of  C6H500,  is  converted  to  C6H50  and  the  dominant  consumption 
channel  of  C6H500  leads  to  C6H402  formation.  The  latter  goes  through  the  sequence  of 
reaction  steps  illustrated  in  Fig.  4.13  that  lead  to  the  formation  of  acetylene. 

In  summary,  a  detailed  chemical  mechanism  for  n-propyl  benzene  oxidation  featuring  1683 
reversible  reactions  and  269  species  has  been  proposed  and  new  reaction  rates  were 
evaluated  and  applied.  The  validation  was  performed  via  comparisons  with  ignition  delay 
times  and  measured  species  concentrations  obtained  by  Eberius  et  al.  [108],  Dagaut  et  al. 
[96]  and  Dagaut  [105]. 


Figure  4.13  Reaction  pathways  for  the  oxidation  of  ^-propyl  benzene  in  a  shock  tube  at  T  =  1448  K, 
$=  1.026  and  P  =  5.06  bar 

The  chemistry  was  partly  developed  based  on  chemical  structure  similarities  with  toluene 
and  propane  for  reactions  occurring  at  the  n-propyl  benzene  branch.  Extensive  comparisons 
between  computed  results  and  the  measurements  for  a  pressure  range  of  1  <  P  (atm)  <10, 
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temperatures  of  900  <  T  (K)  <  1250  and  for  various  stoichiometries  between  0.5  <  O  <  2.0 
illustrate  the  ability  of  the  model  to  predict  the  oxidation  of  n-propyl  benzene. 

Important  reaction  paths  were  identified,  including  the  generation  of  styrene.  The 
significance  of  the  thermal  breakdown  of  the  fuel  at  the  primary  and  secondary  carbon  atom 
was  also  noted  as  small  rate  changes  could  lead  to  large  discrepancies  in  the  intermediate 
species  concentrations.  The  dominant  homolytic  reaction  step  in  both  jet  stirred  reactors  and 
shock  tubes  is  predicted  to  lead  to  benzyl  radical  formation.  The  prediction  and 
measurement  of  large  early  concentrations  of  specific  species  such  as  styrene,  ethyl  benzene 
or  toluene  is  a  clear  indication  that  abstraction  and  homolytic  reactions  occur  at  the 
beginning  of  the  fuel  oxidation  and  promote  the  fuel  decay. 

The  results  presented  here  are  encouraging  and  suggest  that  the  current  reaction  class  based 
approach  can  be  applied  also  to  other  alkyl  substituted  aromatic  fuel  compounds  that  form 
part  of  real  and  surrogate  fuel  formulations.  It  must  also  be  noted  that  as  the  chain  increases 
in  length,  new  categories  of  reaction  steps  are  likely  to  occur  such  as  radical  isomerisations, 
cyclic  transition  states  and  group  shifts  along  the  chain  that  should  be  taken  into  account. 


5.0  The  Oxidation  of  Naphthalene  and  1 -Methyl  Naphthalene 


The  formation  and  oxidation  of  polycyclic  aromatic  hydrocarbons  is  one  of  the  central 
subjects  of  research  in  the  area  of  pollutant  emissions.  Many  studies  in  the  literature  have 
addressed  the  formation  and  growth  of  aromatic  compounds  that  act  as  precursors  to  soot. 
The  control  and  prediction  of  the  combustion  of  aromatics  is  of  direct  importance  to 
combustor  life  times  as  well  as  for  both  environmental  and  health  reasons  as  PAHs  are 
associated  with  tumorigenic  effects.  Hence,  the  ability  to  reproduce  the  chemical  properties 
of  aromatic  hydrocarbon  fuels  is  required. 

The  formation  of  benzene  has  received  particular  attention  and  has  been  subject  to  a 
number  of  studies  [29,  58,  110-113].  Subsequent  molecular  growth  mechanisms, 
specifically  reaction  routes  that  proceed  via  acetylene  based  addition  to  radicals,  have  been 
reported  in  the  literature  [29,  31,  114].  Marinov  et  al.  [42]  and  Castaldi  et  al.  [115]  have 
suggested  that  naphthalene  is  produced  by  the  recombination  of  two  cyclopentadienyl 
radicals  releasing  two  hydrogen  atoms.  The  cyclopentadienyl  recombination  leading  to 
naphthalene  formation  was  also  studied  by  Melius  et  al.  [26]  and  Dean  [116].  Moreover, 
Marinov  et  al.  [117]  and  Richter  et  al.  [118]  highlighted  the  role  of  cyclopentadienyl  radical 
in  premixed  flames.  A  more  recent  study  of  Lindstedt  et  al.  [45]  proposed  a  two-step 
mechanism  involving  the  cyclopentadienyl  route  with  energy  barriers  that  are  in  accordance 
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with  potential  energy  surfaces  determined  by  Melius  et  al.  [26].  The  naphthalene  formation 
and  oxidation  pathways  are  presented  in  this  study  and  compared  with  plug  flow  reactor 
data  obtained  by  Shaddix  [119]. 

5.1  The  Oxidation  of  Naphthalene 

The  reaction  mechanism  derived  for  naphthalene  oxidation  was  based  on  reactions  steps 
and  reaction  classes  resulting  from  studies  of  the  oxidation  of  benzene  and  cyclopentadiene 
with  parameters  adjusted  according  to  molecular  size  and  reactive  site  differences  [54]. 

The  studies  by  Lindstedt  et  al.  [45],  Maurice  [103]  and  Potter  [60]  identified  the  importance 
of  the  linkages  between  aromatic  C5  and  C6  structures  and  the  key  role  of  the 
indene/indenyl  system.  In  particular,  the  oxidation  of  indenyl  was  found  to  exert  a  strong 
influence  on  the  overall  oxidation  kinetics  of  larger  aromatic  species  such  as  naphthalene 
and  1 -methyl  naphthalene  and  on  the  subsequent  product  distribution  of  single  ring 
aromatics.  The  studies  identified  possible  reaction  channels  featuring  molecular  oxygen  and 
provided  preliminary  fits  for  a  limited  temperature  range  (70,  71,  72).  However,  accurate 
estimates  of  the  rates  of  reaction  proved  problematic  [60].  In  the  current  study,  the 
aforementioned  channels  were  considered  along  with  H02  channels  (68,  69)  with  DFT  and 
composite  quantum  mechanical  methods  used  to  calculate  the  potential  energy  surfaces. 
Both  RRKM/ME  theory  and  VTST  approaches  were  subseqently  used  to  derive  estimates 
of  the  rate  constants  [54]. 


c9h7 

+  o2 

=  c6h50 

+  c3h2o 

(71) 

c9h7 

+  o2 

=  c7H7 

+  CO  +  CO 

(70) 

c9h7 

+  o2 

=  c7h60 

+  c2ho 

(72) 

c9h7 

+  ho2 

=  c9h7o 

+  OH 

(68) 

c9h7 

+  ho2 

=  c9h6o 

+  h2 

(69) 

The  full  reaction  mechanism  consists  of  1431  reversible  reactions  and  269  species  and  rate 
parameters  for  all  reactions  discussed  in  the  current  report  are  given  in  Appendix  A. 

The  naphthalene  mechanism  was  validated  under  plug  flow  reactor  conditions  and 
compared  with  measurements  by  Shaddix  [119].  An  offset  of  32  ms  in  the  temporal 
evolution  of  the  experimental  data  was  also  reported  [119].  Predictions  of  reactants  and 
major  species  are  presented  Figs.  5.1  -  5.3.  The  measured  naphthalene  concentrations  show 
scatter  with  the  trend  in  reasonable  agreement  with  computations.  The  fuel  is  decomposed 
via  two  major  channels  that  involve  H  abstraction  via  the  OH  radical  (49%)  and  O  addition 
to  the  ring  leading  to  the  naphthoxy  radical  (44%).  The  rate  assigned  to  reaction  (78)  was 
based  on  the  reaction  rate  of  the  hydrogen  atom  abstraction  from  the  benzene  ring  [58].  The 
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rate  for  reaction  (77)  is  based  on  the  proposal  by  DiNaro  et  al.  [81]  for  the  corresponding 
benzene  reaction. 

C10H8  +  OH  =  C10H7  +  H20  (78) 

C10H8  +  O  =  C10H7O  +  H  (77) 

The  naphthyl  radical  is  further  decomposed  (67%)  to  naphthoxy  via  reaction  (76)  with  the 
assigned  rate  based  on  the  kinetics  of  benzene  as  proposed  by  Frank  et  al.  [95]. 

CioH7  +  02  —  CioH70  +  O  (76) 

The  naphthoxy  radical  decomposes  thermally  via  one  major  pathway  leading  to  the 
formation  of  the  indenyl  radical  and  CO.  The  rate  assigned  to  reaction  (80)  was  based  on 
suggestions  of  Potter  [60]  and  adjusted  from  benzene  kinetics.  The  reaction  is  also  the 
major  (69%)  indenyl  radical  production  channel. 

C10H7O  =  C9H7  +  CO  (80) 

The  indenyl  radical  is  consumed  via  three  major  pathways.  Most  (33%)  recombines  with 
the  hydrogen  atom  to  indene  (65)  with  a  reaction  rate  assigned  based  on  suggestions  of 
Potter  [60].  The  second  route  (19%)  involves  CO  abstraction  via  oxygen  attack  leading  to 
the  formation  of  C8H7  (67).  The  rate  assigned  to  this  channel  was  based  on  suggestions  of 
Maurice  [103].  A  further  17%  of  the  indenyl  concentration  reacts  with  H02  forming  the 
naphthoxy  radical  (68)  with  a  rate  proposed  by  Lindstedt  et  al  [54]. 


c9h7 

+  H 

=  c9H8 

(65) 

c9h7 

+  o 

=  c8H7 

+  CO 

(67) 

c9h7 

+  ho2 

=  c9h7o 

+  OH 

(68) 

Indene  is  formed  via  hydrogen  addition  to  the  indenyl  radical.  Once  it  is  formed,  it 
predominantly  recycles  back  to  the  indenyl  radical.  The  C8H7  radical  decomposes  (100%) 
to  phenyl  acetylene  with  a  rate  adopted  from  Wang  and  Frenklach  [110]. 

C8H7  =  C8H6  +  H  (58) 

Reaction  (58)  is  also  the  major  phenyl  acetylene  formation  channel  (75%).  Phenyl  acetylene 
is  predominantly  oxidised  via  the  reaction  routes  shown  below. 


c8h6 

+  o 

=  C6H5 

+ 

c2ho 

(54) 

c8h6 

+  OH 

=  C8H5 

+ 

h2o 

(52) 

c8h6 

+  o 

=  c8h50 

+ 

H 

(53) 

c8h6 

+  0 

=  c6h50 

+ 

c2h 

(55) 

c8h6 

+  o 

=  c6h5c2o 

+ 

H 

(51) 
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The  major  phenyl  acetylene  consumption  channel  (33%)  occurs  via  oxygen  attack  and 
removes  the  side  chain  via  C-C  cleavage  (54).  The  rate  assigned  to  reaction  (54)  was 
adopted  from  suggestions  of  Lindstedt  et  al.  [45].  The  second  major  consumption  channel 
(52)  occurs  via  H  abstraction  by  OH  attack,  which  is  responsible  for  19%  of  C8H6 
consumption.  The  reaction  was  assigned  a  rate  recommended  by  Frenklach  and  Wang 
[110].  Reactions  (53),  (55)  and  (51)  are  responsible  for  11%,  11%  and  10%  of  the  phenyl 
acetylene  consumption  respectively. 


Figure  5.1  Concentration  profiles  of  major  species  during  naphthalene  oxidation  in  a  plug  flow 
reactor  for  O  =  1.1,  T  =  1197  and  P  =  1  atm.  The  solid  lines  indicate  the  current  computations  and 
the  circles  indicate  the  measurements  [119]. 


Figure  5.2  Concentration  profiles  of  major  species  during  naphthalene  oxidation  in  a  plug  flow 
reactor  for  O  =  1.1,  T  =  1197  and  P  =  1  atm.  The  solid  lines  indicate  the  current  computations  and 
the  circles  indicate  the  measurements  [119]. 

Toluene  is  formed  via  the  benzyl  radical  resulting  from  the  oxidation  of  the  indenyl  radical 
via  molecular  oxygen  attack  (70).  The  rate  assigned  to  reaction  (70)  was  calculated  through 
RRKM/ME  and  VTST  approaches  [54].  Approximately  20%  of  the  benzyl  radical  is 
formed  via  C8H8  oxidation  with  OH  (60)  with  the  rate  adopted  from  Maurice  [103].  A 
further  11%  of  benzyl  radical  is  formed  via  reaction  (34)  with  a  rate  adopted  from  Colket  et 
al  [59]. 
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c9h7 

+  o2 

=  c7H7 

+  CO  +  CO 

(70) 

n 

00 

a 

00 

+  OH 

=  c7H7 

+  CH2o 

(60) 

c7h7 

+  c2h2 

=  c9H8 

+  H 

(34) 

Figure  5.3  Concentration  profiles  of  major  species  during  naphthalene  oxidation  in  a  plug  flow 
reactor  for  O  =  1.1,  T  =  1197  and  P  =  1  atm.  The  solid  lines  indicate  the  current  computations  and 
the  circles  indicate  the  measurements  [119]. 

Benzene  is  formed  (32%)  through  hydrogen  assisted  isomerization  (26)  and  (24%)  from 
benzaldehyde  via  oxygen  addition  and  C02  abstraction  (47).  Benzene  is  a  precursor  to  the 
cyclopentadienyl  radical  from  which  cyclopentadiene  is  formed  (92%).  Benzene  is 
consumed  (47%)  to  the  phenyl  radical  via  H  abstraction  featuring  OH  radical  attack  and 
(40%)  to  the  phenoxy  radical  via  O  addition  and  hydrogen  abstraction.  The  phenyl  radical 
predominantly  (66%)  leads  to  the  phenoxy  radical. 


C6H6(F)  +m 

=  C6H6  +  M 

(26) 

c7h6o  +  0 

-  CgHg  +  C02 

(47) 

The  phenoxy  radical  is  responsible  (77%)  for  cyclopentadienyl  radical  production  via  CO 
expulsion  (28)  with  a  rate  proposed  by  Leung  and  Lindstedt  [58].  The  cyclopentadienyl 
radical  recombines  with  the  H  atom  (92%)  to  the  cyclopentadienyl  radical  (13)  with  a  rate 
adopted  from  Kern  et  al.  [36]. 


c6h5o 

=  C5H5  +  CO 

(28) 

c5h6 

=  C5H5  +  H 

(13) 

Acetylene  is  also  formed  (44%)  via  the  thermal  breakdown  of  C5H5  that  leads  to  C2H2  and 
C3H3.  The  rate  used  for  reaction  (3)  was  proposed  by  Kern  et  al.  [36]  and  was  reduced  by  a 
factor  of  2  due  to  acetylene  overproduction.  The  evolution  of  C2H4  is  also  affected  by  C5H5 
via  reaction  channels  that  include  C2H2  and  C2H20. 

C5H5 


C2H2  +  C3H3 


(3) 
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Vinyl  acetylene  production  is  initially  controlled  by  phenyl  acetylene  consumption  via 
C8H50,  which  thermally  decomposes  to  C7H5  leading  to  C4H4  via  C5H40.  A  rate  by  Ristori 
et  al.  [120]  was  assigned  to  reaction  (15),  which  is  responsible  for  94%  of  the  C4H4 
production. 

C5H40  =  C4H4  +  CO  (15) 

Overall,  the  resulting  mechanism  shows  comparativly  good  agreement  when  evaluated 
against  the  species  profiles  obtained  from  Shaddix  [119]. 

The  work  shows  that  the  naphthalene  oxidation  is  controlled  by  two  major  channels  that 
involve  H  abstraction  via  OH  radical  and  O  addition  to  the  ring.  Most  of  the  oxidation 
process  is  assumed  to  follow  a  similar  behaviour  to  the  benzene  chemistry.  From  the 
comparison  between  the  computations  and  the  measurements,  is  shown  that  the  main 
oxidation  channels  are  well  represented  and  that  the  C9H7  +  02/H02  reactions  have  an 
important  role  in  determining  the  product  distribution  of  single  ring  aromatics. 

5.2  The  Oxidation  of  1 -Methyl  Naphthalene 

Few  studies  on  the  combustion  of  polycyclic  aromatic  fuels  have  been  reported  in  the 
literature  and  in  the  current  work  a  mechanism  for  1 -methyl  naphthalene  was  developed  on 
the  basis  of  reaction  classes  determined  with  toluene  as  a  reference  fuel.  The  choice  is 
natural  since  1 -methyl  naphthalene  is  the  corresponding  bicyclic  analog. 

Shaddix  et  al.  [119,  121,  122]  studied  the  oxidation  of  1 -methyl  naphthalene  in  a  plug  flow 
reactor  and  obtained  species  profiles  under  atmospheric  pressure  conditions  for  a 
temperature  range  of  1165  -  1200  K  and  for  equivalence  ratios  of  0.6,  1.0  and  1.6.  The 
study  arguably  contains  the  first  reported  species  profiles  for  the  oxidation  of  a  polycyclic 
hydrocarbon  fuel.  Based  on  the  obtained  profiles,  a  chemical  mechanism  for  1 -methyl 
naphthalene  was  proposed.  Elevated  pressure  (P  =  13  atm)  ignition  delay  times  were 
measured  by  Pfahl  et  al  [123]  under  shock  tube  conditions  for  a  temperature  range  of  840  - 
1300  K.  Pitsch  [125]  produced  a  reaction  mechanism  capable  of  reproducing  the  ignition 
delay  times  measured  by  Pfahl  et  al.  [123]  as  well  as  ignition  in  the  fuel  lean  mixture 
studied  by  Shaddix  et  al.  [119,  121,  122].  A  more  recent  study  of  the  chemistry  of  1 -methyl 
naphthalene  was  performed  by  Mati  et  al.  [124],  who  studied  the  oxidation  under  jet  stirred 
reactior  conditions  for  a  temperature  range  of  800  <  T  <  1421  K,  a  pressure  range  of  1  <  P  < 
13  atm  and  equivalence  ratios  of  0.5  <  O  <  1.5.  Major  species  concentration  profiles  were 
obtained  by  on-line  GC-MS,  off-line  GC-TCD-FID  and  GC-MS  analyses  and  a  kinetic 
model  was  proposed  and  validated.  Important  reactions  pathways  were  identified  in  the 
above  studies.  Shaddix  et  al.  [122]  highlighted  the  fuel  consumption  paths  that  occur  via  i) 
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abstractions  of  the  benzylic  hydrogen  atom  by  H/O  radicals  or  molecular  oxygen,  ii) 
homolysis  from  C-H  fission  forming  the  1 -naphthyl  methyl  radical,  iii)  displacement  of  the 
methyl  group  by  the  hydrogen  atom  leading  to  the  formation  of  naphthalene  and  iv) 
addition  of  an  oxygen  atom  leading  to  the  formation  of  l-methyl-4-naphthol  or  1 -methyl 
naphthoxy  radical. 

The  oxidation  of  1 -methyl  naphthalene  is  here  studied  and  compared  against  data  from  Mati 
et  al.  [124]  and  Shaddix  [119].  The  starting  point  is  the  detailed  scheme  for  1 -methyl 
naphthalene  derived  by  Potter  [60]  on  the  basis  of  the  chemistry  of  toluene  with 
modifications  and  updates  applied  as  appropriate.  The  reaction  pathways  were  analysed  at 
conditions  corresponding  to  the  JSR  experiments  of  Mati  et  al.  [124]  and  plug  flow  reactor 
conditions  of  Shaddix  [119]  and  Shaddix  et  al.  [121,  122]. 

5.3  The  Oxidation  of  1 -Methyl  Naphthalene  under  JSR  Conditions 

The  1-methyl  naphthalene  mechanism  was  initially  validated  under  jet-stirred  reactor 
conditions  using  species  data  obtained  by  Mati  et  al.  [124].  The  measurements  were  aimed 
at  clarifying  the  decomposition  channels  of  the  fuel.  Three  different  stoichiometries  were 
tested  at  atmospheric  pressure  over  a  temperature  range  of  1090-1450  K  and  at  a  fixed 
residence  time  (x)  of  0.1  s.  The  conditions  are  listed  in  Table  5.1  and  examples  of  species 
profiles  are  shown  in  Figs.  5.4  -  5.15.  The  major  species  detected,  besides  CO  and  C02, 
were  formaldehyde,  methane,  ethane,  ethylene,  acetylene,  cyclopentadiene,  benzene, 
toluene,  styrene,  indene  and  naphthalene. 


P  (atm) 

T  (K) 

X02 

XcilHIO 

0.5 

1.0 

1097-1290 

0.0270 

0.001 

1.0 

1.0 

1094-1400 

0.0135 

0.001 

1.5 

1.0 

1147-1440 

0.0090 

0.001 

Table  5.1  Experimental  and  modelling  conditions  for  the  oxidation  of  1 -methyl  naphthalene  in  a  jet- 
stirred  reactor.  The  species  concentrations  correspond  to  mole  fractions. 


A  reaction  rate  analysis  was  performed  for  a  fuel  lean  mixture  (®  =  0.5)  at  a  temperature  of 
1206  K  at  atmospheric  pressure.  The  computations  show  that  the  major  fuel  decomposition 
pathways  occur  via  hydrogen  abstraction  by  OH  attack  and  oxygen  addition.  The  rates 
assigned  to  the  reactions  were  adopted  from  Potter  [60]  and  determined  on  the  basis  of  the 
toluene  chemistry.  The  OH  attack  on  the  methyl  branch  (86)  is  responsible  for  31%  of  the 
fuel  consumption,  whereas  the  hydrogen  abstraction  from  the  ring  (83)  is  responsible  for 
19%.  The  oxygen  addition  to  the  branch  (88)  accounts  for  12%  of  the  1-methyl  naphthalene 
consumption,  whereas  the  oxygen  addition  to  the  ring  (89)  accounts  for  23%. 
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Figure  5.4  Concentration  profiles  of  intermediate  species  during  1 -methyl  naphthalene  oxidation  in  a 
jet- stirred  reactor  with  O  =  0.5,  P  =  1  atm  and  T  =  1097  -  1290  K.  Circles  are  experimental  data 
[124]  and  the  solid  lines  the  current  simulations. 


Figure  5.5  Concentration  profiles  of  intermediate  species  during  1 -methyl  naphthalene  oxidation  in  a 
jet-stirred  reactor  with  ®  =  0.5,  P  =  1  atm  and  T  =  1097  -  1290  K.  Circles  are  experimental  data 
[124]  and  the  solid  lines  the  current  simulations. 


Figure  5.6  Concentration  profiles  of  intermediate  species  during  1 -methyl  naphthalene  oxidation  in  a 
jet-stirred  reactor  with  ®  =  0.5,  P  =  1  atm  and  T  =  1097  -  1290  K.  Circles  are  experimental  data 
[124]  and  the  solid  lines  the  current  simulations. 
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Figure  5.7  Concentration  profiles  of  intermediate  species  during  1 -methyl  naphthalene  oxidation  in 
jet-stirred  reactor.  ®  =1.0,  P  =  1  atm  and  T  =  1094  -  1400  K.  Circles  are  experimental  data  [124]  and 
the  solid  lines  the  current  simulations. 


Figure  5.8  Concentration  profiles  of  intermediate  species  during  1 -methyl  naphthalene  oxidation  in 
jet-stirred  reactor.  ®  =1.0,  P  =  1  atm  and  T  =  1094  -  1400  K.  Circles  are  experimental  data  [124]  and 
the  solid  lines  the  current  simulations. 


Figure  5.9  Concentration  profiles  of  intermediate  species  during  1 -methyl  naphthalene  oxidation  in 
jet-stirred  reactor.  O  =1.0,  P  =  1  atm  and  T  =  1094  -  1400  K.  Circles  are  experimental  data  [124]  and 
the  solid  lines  the  current  simulations. 
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Figure  5.10  Concentration  profiles  of  intermediate  species  during  1 -methyl  naphthalene  oxidation  in 
jet-stirred  reactor.  ®  =  1.5,  P  =  1  atm  and  T  =  1147  -  1440  K.  Circles  are  experimental  data  [124]  and 
the  solid  lines  the  current  simulations. 


Figure  5.11  Concentration  profiles  of  intermediate  species  during  1 -methyl  naphthalene  oxidation  in 
jet-stirred  reactor.  ®  =  1.5,  P  =  1  atm  and  T  =  1147  -  1440  K.  Circles  are  experimental  data  [124]  and 
the  solid  lines  the  current  simulations. 


Figure  5.12  Concentration  profiles  of  intermediate  species  during  1 -methyl  naphthalene  oxidation  in 
jet-stirred  reactor.  ®  =  1.5,  P  =  1  atm  and  T  =  1147  -  1440  K.  Circles  are  experimental  data  [124]  and 
the  solid  lines  the  current  simulations. 
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CnHio 

+ 

OH 

-  CnH9  + 

h2o 

(86) 

CnHio 

+ 

OH 

=  CnH9P  + 

h2o 

(83) 

CnHio 

+ 

O 

—  ChH90  + 

H 

(88) 

CnHio 

+ 

O 

—  OCnH9  + 

H 

(89) 

The  fuel  is  reformed  from  the  1-methylnaphthyl  radical  (82%)  (-84).  The  reaction  rate  was 
adopted  from  Potter  [60]  and  based  on  the  toluene  kinetics  as  proposed  by  Lindstedt  and 
Maurice  [67].  The  formation  from  HOCnH9  also  makes  a  contribution  (17%)  via  OH 
replacement  (87)  with  HOCnH9  also  formed  from  the  OCnH9  radical  via  H  addition  (92). 
The  OCnH9  radical  is  formed  from  1 -methyl  naphthalene  by  replacement  of  a  hydrogen 
atom  on  the  ring  (89)  with  an  oxygen  atom  (52%).  The  rate  that  stems  from  the  kinetics  of 
toluene  proposed  by  Hoffman  et  al.  [109]  and  from  oxygen  addition  (32%)  to  the  1-methyl- 
4-naphthyl  radical  (CnH9P)  via  molecular  oxygen  attack  (82)  with  a  rate  based  on  the 
toluene  kinetics  proposed  by  Lindstedt  and  Maurice  [67].  The  major  1 -methyl  naphthalene 
decomposition  and  formation  pathways  are  shown  in  Figs.  5.13  -  5.14. 


CnHio 

=  c„H9 

+  H 

(84) 

CnHio 

+ 

OH 

=  HOCnH9 

+  H 

(87) 

OCnH9 

+ 

H 

=  HOCnH9 

(92) 

CnHio 

+ 

0 

=  OCuH9 

+  H 

(89) 

ChH9P 

+ 

o2 

=  OCuH9 

+  0 

(82) 

Naphthalene  is  formed  via  thermal  decomposition  of  the  CnH90  radical  (91)  (65%),  which 
is  formed  from  oxygen  attack  on  1-methyl  naphthalene  (88)  (Fig.  5.15). 

CnH10  +0  =  CnH90  +  H  (88) 

CnH90  =  C10H8  +  CHO  (91) 

Once  naphthalene  is  formed,  hydrogen  abstraction  via  OH  radical  attack  (78),  leading  to  the 
formation  of  the  naphthyl  radical  contributes  49%  to  consumption  of  the  former.  Moreover, 
43%  of  the  naphthalene  forms  naphthaldehyde  (77)  by  oxygen  addition  and  hydrogen 
abstraction.  The  naphthalene  decomposition  reactions  were  assigned  rates  based  on  the 
corresponding  benzene  channels.  The  naphthaldeyde  thermally  decomposes  by  CO 
abstraction  leading  to  the  formation  of  the  indenyl  radical  and  indene. 

C10H8  +  O  =  C10H7O  +  H  (77) 

C10H8  +  OH  =  C10H7  +  H20  (78) 

The  formation  of  indene  occurs  (98%)  from  the  indenyl  radical.  The  indenyl  radical  is 
formed  in  a  sequence  of  reaction  steps  initiated  by  the  formation  of  the  OCnH9  radical.  The 
latter  thermally  decomposes  via  CO  expulsion  from  the  ring  leading  to  the  formation  of 
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C9H6CH3.  The  latter  species  recombines  with  a  hydrogen  atom  leading  to  the  formation  of 
C9H7CH3,  which  thermally  decomposes  to  indenyl  and  methyl  (50%).  The  rates  of  the  last 
two  steps  were  assigned  from  suggestions  made  by  Laskin  and  Lifshitz  [126].  In  addition, 
the  indenyl  radical  is  formed  by  the  thermal  decomposition  of  naphthaldehyde  (16%).  A 
more  direct  link  to  the  fuel  is  provided  by  (12%)  the  indenyl  radical  and  acetylene 
formation  from  the  1-methylnaphthyl  radical  (81).  The  reverse  reaction  step  was  also  tested 
with  a  rate  relative  to  the  benzyl  radical  formation  channel  proposed  by  Colket  and  Seery 
[62],  but  found  to  lead  to  excessive  acetylene  and  indene  concentrations.  By  contrast,  the 
use  of  the  1 -methyl  naphthtyl  radical  thermal  decomposition  channel  (81)  with  a  rate 
relative  to  the  benzyl  radical  thermal  decomposition  proposed  by  Braun-Unkhoff  [77]  leads 
to  excellent  agreement  between  computed  and  measured  concentration  profiles  for  both 
acetylene  and  indene. 

CnH9  =  C9H7  +  C2H2  (81) 

The  oxidation  of  the  indenyl  radical  by  H02  leads  to  C9H70  formation  (68).  Density 
Functional  Theory  and  composite  quantum  mechanical  methods  were  used  to  calculate  the 
potential  energy  surfaces  for  the  channel.  Both  RRKM/ME  theory  and  VTST  were 
subsequently  used  to  derive  an  estimate  of  the  rate  constant  [54].  The  species  is  a  precursor 
to  phenyl  acetylene  (C8H6)  formed  through  C8H7  decomposition.  A  rate  proposed  by  Pitsch 
[125]  was  adopted  for  the  CO  expulsion  from  C9H70  (75).  The  thermal  decomposition  of 
C8H7  is  responsible  for  71%  of  the  phenyl  acetylene  formation  (58)  with  a  rate  based  on 
suggestions  of  Wang  and  Frenklach  [110]. 


c9h7  +  H02  = 

c9h7o 

+ 

OH 

(68) 

c9h7o 

c8h7 

+ 

CO 

(75) 

c8h7 

c8h6 

+ 

H 

(58) 

Computations  show  that  molecular  oxygen  attack  (70)  forms  another  important  indenyl 
oxidation  step  responsible  for  50%  of  the  benzyl  radical  formation  with  the  rate  determined 
by  Lindstedt  et  al.  [54].  The  benzyl  radical  leads  to  toluene  via  reaction  (38)  with  the  rate 
determined  by  Oehlschlaeger  et  al.  [80]. 

C9H7  +  02  C7H7  +  CO  +  CO  (70) 

C7H8  =  C7H7  +  H  (38) 

In  addition  to  the  fulvene  reaction  sequence,  25%  of  benzene  is  formed  from  benzaldehyde 
and  a  rate  proposed  by  Potter  [60]  was  used  for  reaction  (47). 

-  CgHg  +  C02 


c7h6o  +  O 


(47) 
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Benzene  is  mostly  formed  (35%)  from  fulvene  via  isomerization  (26)  assigned  a  rate 
proposed  by  Marinov  et  al.  [63].  Fulvene  is  produced  (100%)  by  the  methyl  cyclo 
pentadienyl  radical  after  thermal  decomposition  (27)  and  assigned  a  rate  proposed  by 
Lindstedt  and  Rizos  [46]. 


C5H4CH3  =  C6H6(F)  +  H 

C6H6(F)  +  M  =  C6H6  +  M 


(27) 

(26) 


C11H9O 


Figure  5.1  Major  1 -methyl  naphthalene  decomposition  pathways. 

Benzene  is  consumed  by  reactions  leading  to  phenyl  (42%)  and  phenoxy  (48%)  radicals. 
The  phenyl  radical  also  predominantly  (65%)  leads  to  phenoxy  which  in  turn  leads  to  the 
formation  of  the  cyclopentadienyl  radical  (80%)  by  CO  expulsion.  The  cyclopentadienyl 
radical  follows  two  major  decomposition  routes;  28%  thermally  decomposes  to  acetylene 
and  the  propargyl  radical  and  25%  recombines  with  hydrogen  leading  to  cyclopentadiene. 
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Cyclopentadiene  is  predominantly  (90%)  produced  by  the  recombination  of  the 
cyclopentadienyl  radical  with  atomic  hydrogen.  A  rate  proposed  by  Kern  et  al.  [36]  was 
applied  to  this  channel.  Cyclopentadiene  then  recycles  back  to  cyclopentadienyl  radical. 

C5H6  =  C5H5  +  H  (13) 

Acetylene  is  formed  via  two  major  thermal  decomposition  channels.  The  dominant  (39%)  is 
the  thermal  decomposition  of  the  1 -methyl  naphthyl  radical  that  forms  acetylene  and  the 
indenyl  radical  (81)  and  acetylene  is  also  formed  (21%)  via  the  thermal  breakdown  of 
cyclopentadienyl  radical  (3). 

CnH9  =  C9H7  +C2H2  (81) 

C5H5  =  C3H3  +  C2H2  (3) 

The  rate  of  Kern  et  al.  [36]  was  applied  to  reaction  channel  (3)  but  found  to  produce 
excessive  acetylene  concentrations,  hence  an  adjustment  was  made  by  reducing  the  rate  by 
50%. 


HOC11H9  OC11H9 


Figure  5.2  Major  1 -methyl  naphthalene  formation  pathways. 
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CnHio  ChH90  CioHs 

Figure  5.3  Naphthalene  formation  pathways 


5.4  The  Oxidation  of  1 -Methyl  Naphthalene  under  TFR  Conditions 

Shaddix  [119]  performed  gas-phase  sampling  to  study  the  oxidation  of  1 -methyl 
naphthalene  in  the  Princeton  turbulent  flow  reactor  and  obtained  concentration  profiles  for 
major  species  at  atmospheric  pressure  and  at  a  temperature  of  approximately  1 170  K.  Major 
fuel  consumption  routes  and  reaction  classes  were  identified.  Computations  were  performed 
corresponding  to  the  experimental  conditions  for  stoichiometric  and  fuel  rich  mixtures  as 
shown  in  Table  5.2. 


<D 

T  (K) 

P  (atm) 

[02]  (ppm) 

[CnH10]  (ppm) 

1.0 

1169 

1.0 

14850 

1100 

1.5 

1166 

1.0 

9900 

1100 

1.5 

1198 

1.0 

9900 

1100 

Table  5.2  Experimental  and  modelling  conditions  for  the  oxidation  of  1 -methyl  napthalene  in  a 
turbulent  flow  reactor. 


Selected  comparisons  between  computations  and  experimental  data  are  shown  in  Figs.  5.16 
-5.19.  The  fuel  profile  is  again  well  reproduced,  suggesting  that  the  1 -methyl  naphthalene 
sub-mechanism  is  adequate.  Moreover,  naphthalene  and  indene,  which  are  two  of  the  major 
initial  products  of  the  fuel  decomposition  process,  are  very  well  reproduced  for  all  the  three 
test  cases  leading  to  good  agreement  between  computations  and  measurements  for 
subsequent  products  such  as  phenylacetylene,  styrene  and  toluene. 

The  direct  removal  of  the  methyl  group,  either  by  replacement  or  via  thermal 
decomposition,  is  a  minor  process  contributing  <  5%  to  the  fuel  breakdown.  Once  the  1- 
methylnaphthyl  radical  is  formed,  the  proposed  thermal  decomposition  channel  leading  to 
the  indenyl  radical  and  acetylene  is  important  and  affects  the  growth  and  evolution  of  major 
intermediate  species  such  as  indene,  styrene,  phenylacetylene,  toluene,  benzene, 
cyclopentadiene  and  acetylene. 
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In  summary,  the  level  of  agreement  obtained  between  the  computed  concentrations  and 
measurements  the  naphthalene  and  indene  implies  that  the  accuracy  of  the  fuel  decay  paths 
that  involve  oxygen  addition  and  lead  to  the  formation  of  CnH90  and  OCnH9is  adequate. 
The  further  decomposition  of  the  oxygenated  species  leads  to  the  formation  of  Ci0H8  and 
significant  amounts  of  C9H7.  The  branching  ratios  between  abstraction  and  addition 
pathways  are  also  important  for  1 -methyl  naphthalene  combustion  at  intermediate 
temperatures  and,  again,  the  level  of  agreement  obtained  is  encouraging.  Furtermore,  the 
distribution  of  the  single-ring  aromatics  is  also  well  reproduced. 


Figure  5.16  Concentration  profiles  of  major  species  during  1 -methyl  naphthalene  oxidation  in  a  flow 
reactor  (®  =  1.0,  P  =  1  atm,  T  =  1169  K).  Circles  indicate  measurements  [119]  and  solid  lines  the 
current  computations. 


Figure  5.17  Concentration  profiles  of  major  species  during  1 -methyl  naphthalene  oxidation  in  a  flow 
reactor  (®  =  1.0,  P  =  1  atm,  T  =  1169  K).  Circles  indicate  measurements  [119]  and  solid  lines  the 
current  computations. 
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Figure  5.18  Concentration  profiles  of  major  species  during  1 -methyl  naphthalene  oxidation  in  a  flow 
reactor  (®  =  1.0,  P  =  1  atm,  T  =  1169  K).  Circles  indicate  measurements  [119]  and  solid  lines  the 
current  computations 


Figure  5.19  Concentration  profiles  of  major  species  during  1 -methyl  naphthalene  oxidation  in  a  flow 
reactor  (®  =  1.0,  P  =  1  atm,  T  =  1169  K).  Circles  indicate  measurements  [119]  and  solid  lines  the 
current  computations. 


6.0  Conclusions 


The  main  objective  of  the  present  work  has  been  the  development  and  validation  of  gas 
phase  chemical  kinetic  mechanisms  for  aromatic  components  of  aviation  fuels  and 
surrogates.  A  reaction  class  based  approach  was  used  for  the  generation  of  the  chemical 
mechanisms  and  a  broad  range  of  reactors  and  conditions  were  considered.  The 
mechanisms  feature  molecules  up  to  Cn  that  involve  aliphatic  compounds  and  sub¬ 
mechanisms  for  mono-substituted  and  polyaromatic  hydrocarbons.  The  starting  point  for 
the  work  was  the  development  and  validation  of  a  mechanism  for  cyclopentadiene,  which 
constitutes  a  major  intermediate  component  in  the  combustion  of  higher  aromatics.  The 
work  was  subsequently  extended  to  toluene  and  n-propyl  benzene  and  two-ringed  aromatics 
such  as  naphthalene  and  1 -methyl  naphthalene. 
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The  oxidation  of  cyclopentadiene  was  studied  under  plug  flow  reactor  conditions  and 
concentration  profiles  were  computed  and  compared  with  measurements  of  Butler  [55]. 
Reaction  rates  for  channels  featuring  O,  H02  and  OH  were  updated  and  validated.  A  variety 
of  cases  varying  from  stoichiometric  to  fuel  rich  were  tested  and  reaction  rate  analyses 
performed  for  representative  cases.  The  pyrolysis  of  cyclopentadiene  was  also  studied  for  a 
variety  of  fuel  concentrations.  The  detection  of  high  concentrations  of  species  such  as 
naphthalene,  toluene,  benzene  and  acetylene  suggest  pathways  leading  to  poly-aromatic 
hydrocarbons. 

The  combustion  of  toluene  was  further  analysed  under  shock  tube  conditions.  Pyrolytic  and 
oxidation  conditions  corresponding  to  those  obtained  by  Braun-Unkhoff  et  al.  [77], 
Vasudevan  et  al.  [70]  and  Burcat  et  al.  [69]  were  included.  Reaction  rates  were  updated  and 
new  reaction  pathways  added  to  the  starting  mechanism.  Hydrogen  atom  profiles  were 
calculated  under  pyrolytic  conditions  and  shown  to  capture  the  measured  profiles  of  Braun- 
Unkhoff  et  al.  [77]  reasonably  well.  The  implication  is  that  the  applied  branching  ratio  of 
the  toluene  thermal  decomposition  to  the  benzyl  radical  or  phenyl  and  methyl  radicals  is 
arguably  correct. 

Time  resolved  OH  profile  was  also  computed  under  shock  tube  conditions  and  compared 
against  experimental  data  of  Vasudevan  et  al.  [70].  Reaction  rate  analyses  were  performed 
and  important  pathways  identified.  The  results  highlight  the  toluene  thermal  decomposition 
to  the  benzyl  radical  and  phenyl  and  methyl  radicals,  as  well  as  hydrogen  assisted  reactions 
leading  to  the  production  of  the  benzyl  radical  and  benzene.  The  results  suggest  that  the 
reaction  class  based  approach  can  be  applied  to  other  methyl-substituted  aromatic  fuel 
components  that  form  part  of  real  and  surrogate  fuel  formulations. 

A  reaction  mechanism  for  n-propyl  benzene  was  developed  on  the  basis  of  the  behaviour  of 
benzene,  toluene  and  propane.  The  mechanism  was  validated  under  jet-stirred  reactor 
conditions  for  both  low  and  high  pressures  for  a  wide  range  of  stoichiometries  and  against 
ignition  delay  times.  A  reaction  rate  analysis  highlighted  the  importance  of  the  hydrogen 
abstraction  reactions,  which  generate  styrenes  and  promote  the  fuel  consumption.  The 
thermal  breakdown  of  the  n-propyl  benzene  chain  at  the  primary  and  secondary  sites  is 
significant  and  small  rate  changes  can  cause  big  discrepancies  in  the  intermediate  species 
concentrations.  Large  early  concentrations  of  styrene,  ethyl  benzene  and  toluene  provide  a 
clear  indication  that  homolytic  reactions  occur  and  promote  the  fuel  decay. 

The  work  was  expanded  to  two-ring  aromatic  compounds  such  as  naphthalene  and  tested 
under  plug  flow  reactor  conditions  at  atmospheric  pressure  and  compared  to  measurements 
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obtained  by  Shaddix  [119].  Reaction  rates  for  the  C9H7  +  02/H02  were  obtained  by 
Lindstedt  et  al.  [54]  and  a  reaction  rate  analysis  performed  showed  that  the  channels  play  a 
significant  role  to  the  product  distribution  of  single  ring  aromatics.  It  was  also  shown  that 
hydrogen  abstraction  reaction  via  OH  and  oxygen  addition  to  the  ring  constitute  major  fuel 
consumption  pathways.  However,  due  to  limited  availability  of  experimental  data  in  the 
literature,  only  a  single  test  case  was  analysed  in  this  work.  Hence,  measurements  of  wider 
range  of  conditions  are  desirable  for  further  refinements.  However,  the  much  more 
extensive  work  performed  on  the  1 -methyl  naphthalene  system,  outlined  below,  is 
complementary  and  adds  further  confidence  in  the  modelling. 

The  1 -methyl  naphthalene  sub-mechanism  was  based  on  the  toluene  chemistry  due  to 
structural  similarities.  Reaction  rates  were  updated  and  reaction  steps  added  to  the 
mechanism.  Validation  was  performed  under  plug  flow  and  jet  stirred  reactor  conditions  at 
atmospheric  pressure  for  a  variety  of  stoichiometries.  Results  were  compared  against 
measurements  obtained  by  Shaddix  and  co-workers  [119,  121]  and  Mati  et  al.  [124].  A 
reaction  rate  analysis  was  performed  for  both  plug  flow  and  jet  stirred  reactor  conditions 
and  showed  that  the  fuel  decomposes  via  OH  hydrogen  abstraction  from  the  chain  and  the 
ring  and  oxygen  addition  to  both  the  chain  and  the  ring.  However,  the  attack  on  the  chain 
side  dominates  the  consumption  of  the  fuel.  The  proposed  thermal  decomposition  of  the  1- 
methylnaphthyl  radical  to  the  indenyl  radical  and  acetylene,  was  shown  to  constitute  a 
major  channel  for  the  evolution  of  intermediate  species.  Moreover,  the  reaction  sequence  of 
C9H7  +  02/H02  was  again  shown  to  be  important  for  the  distribution  of  single-ring 
aromatics.  The  good  agreement  between  the  computations  and  the  measurements  for  the 
current  1 -methyl  naphthalene  is  encouraging. 

Overall,  the  current  work  has  resulted  in  a  significantly  improved  understanding  of  key 
reaction  pathways  for  both  the  oxidation  and  growth  of  aromatic  species  in  combustion 
applications  and  opens  up  a  route  for  further  improvements  in  the  predictions  of  not  only 
global  properties,  such  as  ignition  delay  times,  but  also  the  speciation  characteristics  that 
are  essential  for  the  prediction  of  pollutant  emissions. 
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Appendix  A 

3.3  Reaction  Mechanisms 

The  detailed  chemical  reaction  mechanism  that  was  developed  and  utilized  in  the 
present  work  is  shown  in  Table  A.l.  Reaction  rates  are  expressed  in  the  form: 


k  =  ATnexp(-Ea/RT)  (A.l) 

where  A  =  frequency  factor,  (kmol/m3)(l  m)/s,  m  is  the  order  of  the  reaction 
n  =  temperature  dependence  exponent 
Ea  =  activation  energy,  kJ/mole 
R  =  ideal  gas  constant,  8.314  kJ/K/kmol 
T  =  temperature,  K 

Unimolecular  and  recombination  reactions  generally  exhibit  complex  pressure  dependence. 
Troe  [12]  proposed  a  generalised  expression  for  the  fall-off  behaviour  of  such  reactions: 


ka\M\  +  /<., 


(A.2) 


where  k^  is  the  high  pressure  limit  rate  constant,  k0  is  the  low  pressure  limit,  [M]  is  the 
concentration  of  the  third  body  collision  partner  and  F  is  the  broadening  factor  which  is 
defined  as  : 


logF  = 


_ log^c _ 

\  +  {og(k0[M]/k„)/N}2 


(A. 3) 


where  N  =  0.75  -  1.271ogFc 


(A. 4) 


and  Fc  is  a  temperature  dependent  parameter  specifically  defined  for  a  particular  reaction. 
The  reaction  mechanism  follows. 
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TABLE  A.1 

Rate  parametyers  for  elementary  reactions  discussed  in  the  current  report 
(units  are  kmol,  m3,  s,  K  and  kJ/mol) 


No 

Reaction 

A 

n 

Ea 

Ref 

1 

C5H5  + 

ch3 

=  C5H4CH3+  H 

2.0000E+10 

0.00 

26.00 

[46] 

2 

C5H5  + 

C2H2 

=  C7H7 

3.7200E+08 

0.00 

34.72 

[81] 

3 

C5H5 

=  C3H3  +  C2H2 

3.1500E+13 

0.075 

260.66 

Adj.[36] 

4 

C5H5  + 

0 

=  C4H5(T)+  CO 

2.1260E+53 

-12.56 

89.96 

[15] 

5 

C5H5  + 

HC>2 

=  C5H5O  +  OH 

6.2480E+62 

-16.16 

82.82 

[15] 

6 

C5H5  + 

C5H5 

=  C10H9F  +  H 

2.0000E+10 

0.00 

34.00 

[45] 

7 

C5H5  + 

H02 

=  C5H4O  +  H20 

6.2480E+62 

-16.16 

82.82 

[15] 

8 

C5H5 

=  c5h5(L) 

1.0740E+38 

-7.49 

307.55 

[15] 

9 

c5h6  + 

H 

=  C5H5  +  h2 

3.4970E+04 

18.56 

18.74 

[15] 

10 

c5h6  + 

0 

=  C5H5  +  OH 

1.8100E+10 

0.00 

12.87 

[58] 

11 

C5H6  + 

OH 

=  c5h5  +  H20 

1.1430E+06 

1.18 

-1.87 

[57] 

12 

C5H6  + 

C5H5 

=  CgHs  +  CH3 

3.0000E+13 

0.00 

177.10 

[25] 

13 

c5h6 

=  C5H5  +  H  15 

1 .0000E+19 

-0.65 

368.19 

[36] 

14 

C5H7 

=  C5H6  +  H 

3.1600E+15 

0.00 

185.40 

[46] 

15 

C5H4O 

=  C4H4  +  CO 

2.5000E+1 1 

0.00 

221.90 

Adj.  [120] 

16 

c6h5  + 

c2h4 

=  C8H8  +  H 

2.5100E+09 

0.00 

25.93 

[193] 

17 

c6h5  + 

H 

=  c6h4  +  h2 

1 .5000E+1 1 

0.00 

0.00 

[58] 

18 

c6h5  + 

H 

=  C6Hs 

7.8300E+10 

0.00 

0.00 

[91] 

19 

c6h5  + 

O2 

=  C6H50  +  0 

2.6000E+10 

0.00 

25.61 

[95] 

20 

c6h5  + 

02 

=  C6H500 

3.0800E+09 

-0.15 

0.66 

[60] 

21 

C6H6  + 

C2H3 

=  C8H8  +  H 

7.4900E+08 

0.00 

26.77 

[45] 

22 

C6H6  + 

H 

=  C6H5  +  h2 

2.5000E+1 1 

0.00 

66.94 

[58] 

23 

C6H6  + 

0 

=  C6H50  +  H 

2.4000E+10 

0.00 

19.53 

[81] 

24 

+ 

CO 

X 

CO 

0 

0 

X 

O 

+ 

LO 

X 

CO 

0 

II 

2.0000E+10 

0.00 

61.52 

[58] 

25 

+ 

CO 

X 

CO 

0 

OH 

=  c6h5  +  H20 

1.6300E+05 

1.42 

6.10 

[58] 

26 

C6H6(F)  h 

■  M 

=  C6H6  +  M 

3.0000E+09 

0.500 

8.37 

[63] 

27 

C5H4CH3 

=  C6H6(F)+  h 

1.0000E+14 

0.00 

217.00 

[46] 

28 

C6H50 

=  C5H5  +  CO 

4.5000E+1 1 

0.00 

126.68 

[58] 

29 

C6H50  h 

-  H 

=  CeHsOH 

2.5000E+1 1 

0.00 

0.00 

[81] 

30 

C6H4O2 

+  0 

=  C6H302  +  OH 

1.4000E+10 

0.00 

61.55 

[81] 

31 

C6H500 

=  C6H50  +  0 

4.2700E+15 

-0.70 

138.27 

[81] 
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No 

Reaction 

A 

n 

Ea 

Ref 

32 

C6H500 

=  C6H402  +  H 

4.0000E+08 

0.00 

0.00 

[81] 

33 

C7H7 

+ 

0 

=  C7H60  +  H 

3.5000E+10 

0.00 

0.00 

[197] 

34 

C7H7 

+ 

C2H2 

=  C9H8  +  H 

1.0000E+09 

0.00 

20.92 

[59] 

35 

C7H7 

+ 

02 

=  C7H7O  +  0 

8.6400E+10 

0.00 

131.36 

[60] 

36 

C7H7 

+ 

02 

=  C7H7OO6  k00  = 

1.3400E+06 

1.10 

0.00 

[60] 

k0  = 

5.3100E+19 

-3.30 

0.00 

37 

C7H7 

X 

+ 

CD 

X 

0 

II 

8.2000E+14 

0.00 

337.52 

[80] 

38 

C7H8 

=  C7H7  +  H 

2.0900E+15 

0.00 

366.14 

[80] 

39 

C7H8 

=  C6H5  +  CH3 

2.6600E+16 

0.00 

409.52 

[80] 

40 

C7H8 

+ 

H 

=  C7H7  +  h2 

1.2600E+12 

0.00 

61.99 

[80] 

41 

c7h8 

+ 

H 

=  C6H6  +  CHs 

5.7800E+10 

0.00 

33.84 

[80] 

42 

c7h8 

+ 

H 

=  C7H7P  +  h2 

2.5000E+1 1 

0.00 

66.94 

[57] 

43 

c7h8 

+ 

OH 

=  C7H7  +  H20 

5.1900E+06 

1.00 

3.65 

[91] 

44 

C7H8 

+ 

CHs 

=  CH4  +  C7H7 

3.1600E+09 

0.00 

0.00 

[80] 

45 

C7H8 

+ 

0 

=  OC7H7  +  H 

3.1000E+10 

0.00 

16.62 

[109] 

46 

C7H8 

+ 

0 

=  C7H7  +  OH 

6.3000E+08 

0.00 

0.00 

[109] 

47 

c7h6o 

+ 

0 

=  CeH6  +  CO2 

2.0000E+10 

0.00 

0.00 

[60] 

48 

C7H70H  + 

c6h5 

=  C7H60  +  C6H6  +  H 

1.4000E+09 

0.00 

18.41 

[39] 

49 

C7H700 

=  C7H60  +  OH 

1.0000E+10 

0.00 

121.00 

[103] 

50 

C7H700 

=  C7H7O  +  0 

7.8300E+16 

0.00 

244.78 

Est.(31) 

51 

c8h6 

+ 

0 

=  c6h5c2o+  h 

2.1900E+03 

2.09 

6.53 

[60] 

52 

c8h6 

+ 

OH 

=  C8H5  +  H20 

2.1000E+10 

0.00 

19.10 

[61] 

53 

c8h6 

+ 

0 

=  C8H50  +  H 

2.4000E+10 

0.00 

19.53 

[60] 

54 

c8h6 

+ 

0 

=  C6H5  +  C2HO 

6.5100E+03 

2.09 

6.54 

[45] 

55 

c8h6 

+ 

0 

=  c6h5o  +  c2h 

2.2000E+10 

0.00 

18.95 

est. 

56 

c8h7 

+ 

OH 

=  C6H5  +  CHsCO 

3.0000E+10 

0.00 

0.00 

[103] 

57 

c8h7 

+ 

0 

=  C6H5  +  C2H2O 

3.0000E+10 

0.00 

0.00 

[103] 

58 

c8h7 

=  C8H6  +  H 

1.3000E+41 

-8.65 

46.10 

[110] 

59 

c8h8 

+ 

OH 

=  C8H7(P)  +  H20 

2.1000E+10 

0.00 

19.10 

[103] 

60 

c8h8 

+ 

OH 

=  C7H7  +  CH20 

3.0000E+10 

0.00 

0.00 

[103] 

61 

c8h8 

+ 

0 

=  C7H7  +  CHO 

4.2020E+03 

1.88 

0.76 

[57] 

62 

c8h8 

+ 

0 

=  C6H5  +  CHsCO 

2.5000E+03 

1.88 

0.76 

[60] 

63 

C8Hio 

+ 

H 

=  c8h9  +  h2 

1.2600E-01 

3.44 

13.05 

[103] 

64 

C8Hio 

+ 

H 

=  CeH6  +  C2H5 

1.2000E+10 

0.00 

21.33 

[107] 

65 

c9h7 

+ 

H 

=  C9H8 

2.0000E+1 1 

0.00 

0.00 

[60] 

66 

c9h7 

+ 

0 

=  C8H7(P)+  CO 

4.5000E+10 

0.00 

0.00 

[60] 
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No 

Reaction 

A 

n 

Ea 

Ref 

67 

C9H7  +  0 

O 

O 

+ 

X 

co 

O 

II 

1.0000E+11 

0.00 

0.00 

[103] 

68 

C9H7  +  HO2 

=  C9H7O  +  OH 

8.5770E+65 

-16.69 

101.78 

[54] 

69 

C9H7  +  HO2 

=  C9H60  +  H20 

1.9600E+30 

-6.14 

228.50 

[54] 

70 

C9H7  +  O2 

=  C7H7  +  CO+  CO 

1.7420E+09 

0.31 

123.10 

[54] 

71 

C9H7  +  O2 

=  C6H50  +  C3H2O 

1.8760E+10 

-0.05 

124.68 

[54] 

72 

C9H7  +  O2 

=  C7H60  +  C2HO 

1.7420E+09 

0.31 

123.10 

[54] 

73 

C9H8  +  H 

=  C9H7  +  h2 

2.1 120E+10 

0.00 

9.45 

[60] 

74 

C9H8  +  OH 

=  C9H7  +  H20 

8.6230E+05 

1.18 

-1.87 

[60] 

75 

C9H7O 

=  C8H7  +  CO 

2.0000E+14 

0.00 

155.00 

[125] 

76 

C10H7  +  O2 

=  C10H7O  +  0 

2.1500E+10 

0.00 

25.61 

[60] 

77 

CioH8  +  0 

=  C10H7O  +  H 

2.5000E+10 

0.00 

19.54 

[60] 

78 

CioH8  +  OH 

=  C10H7  +  H2O 

1.7000E+05 

1.42 

6.07 

[60] 

79 

C10H9F 

=  CioH8  +  H 

3.0000E+13 

0.00 

197.00 

[103] 

80 

C10H7O 

11 

O 

CO 

X 

+ 

0 

0 

1.8000E+11 

0.00 

183.68 

[60] 

81 

C11H9 

=  C9H7  +  C2H2 

8.0000E+09 

0.00 

186.60 

Adj.[77] 

82 

C11H9P  +  O2 

II 

0 

p 

X 

CO 

+ 

0 

1.9575E+10 

0.00 

25.61 

[60] 

83 

C11H10  +  OH 

=  C11H9P  +  H20 

1.7000E+05 

1.42 

6.07 

[60] 

84 

C11H10 

=  C11H9  +  H 

4.4800E+15 

0.00 

381.30 

[60] 

85 

CM 

O 

+ 

O 

X 

6 

=  C11H9  +  H02 

1.2670E+09 

0.00 

166.28 

[60] 

86 

C11H10  +  OH 

=  C11H9  +  H20 

4.1520E+06 

1.00 

3.66 

[60] 

87 

C11H10  +  OH 

=  HOC11H9  +  H 

1.8320E+09 

0.00 

-1.49 

[60] 

88 

C11H10  +  0 

=  C11H9O  +  H 

1.2420E+10 

0.00 

16.62 

[60] 

89 

C11H10  +  0 

=  OC11H9  +  H 

2.4000E+10 

0.00 

16.54 

[60] 

90 

ChH80  +  0 

=  CioH8  +  CO2 

1.6000E+10 

0.00 

0.00 

[60] 

91 

C11H9O 

=  CioH8  +  CHO 

4.0000E+12 

0.00 

79.00 

[60] 

92 

OC11H9  +  H 

=  HOC11H9 

2.0000E+1 1 

0.00 

0.00 

[60] 
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TABLE  A.2 

Sub-mechanism  for  n-propyl  benzene  oxidation 
(units  are  kmol,  m3,  s,  K  and  kJ/mol) 


No 

Reaction 

A 

n 

Ea 

Ref 

93 

^gHio 

+  H 

= 

C6H5  +  C3H6 

1.0000E+10 

0.00 

44.57 

[211] 

94 

’CsHio 

+  H 

= 

c8h8  +  ch3 

3.9800E+10 

0.00 

9.97 

[211] 

95 

’CsHio 

+  H 

= 

c6h6  +  C3H5(S) 

1.5848E+11 

0.00 

52.35 

[211] 

96 

’CsHio 

+  H 

= 

C9H9(T)+  h2 

1.5848E+11 

0.00 

46.53 

[211] 

97 

^gHio 

+  0 

= 

c7h6o  +  C2H4 

5.2300E+04 

1.57 

-2.63 

[60] 

98 

^gHio 

+  02 

= 

C9H9CO  +  H02 

2.0000E+10 

0.00 

184.10 

[60] 

99 

^gHio 

+  02 

= 

C7H60  +  CH3CHO 

1.0000E+10 

0.00 

0.00 

[60] 

100 

’CsHio 

+  HC>2 

= 

C7H60  +  C2H4  +  OH 

1.0000E+09 

0.00 

0.00 

[60] 

101 

’CsHio 

+  OH 

= 

C7H60  +  C2H5 

1.0000E+10 

0.00 

0.00 

[60] 

102 

2C9Hio 

+  H 

= 

3C9Hii 

7.2300E+09 

0.00 

12.14 

[60] 

103 

2CgHio 

+  H 

= 

C6H5  +  C3H6 

1.0000E+11 

0.00 

44.57 

[60] 

104 

o 

X 

O) 

o 

CM 

+  H 

= 

C7H7  +  C2H4 

2.6000E+05 

1.50 

8.36 

Present 

Work 

105 

2CgHio 

+  H 

= 

C6Hs  +  C3H5(S) 

1.5848E+11 

0.00 

52.35 

[60] 

106 

o 

X 

O) 

o 

CM 

+  O 

= 

c8h8  +  CH2O 

5.2300E+04 

1.57 

2.63 

[60] 

107 

'CsHu 

= 

’CsHio  +  H 

6.3000E+13 

0.00 

154.39 

[60] 

108 

'CsHu 

= 

C7H7  +  C2H4 

2.0000E+10 

0.00 

123.50 

[60] 

109 

2C9Hh 

= 

2C9H10  +  H 

3.1500E+13 

0.00 

154.39 

[60] 

110 

2C9Hh 

= 

C8H8  +  ch3 

2.0000E+10 

0.00 

123.43 

[60] 

111 

2C9Hh 

= 

^gHio  +  H 

3.1500E+13 

0.00 

154.39 

[60] 

112 

2C9Hh 

= 

C6H6  +  C3H5(S) 

4.0000E+13 

0.00 

197.75 

[60] 

113 

2C9Hh 

= 

C6H6  +  C3H5(A) 

0.00 

0.00 

0.00 

[60] 

114 

3C9Hh 

+  H02 

= 

C8Hio  +  CHO  +  OH 

2.5000E+09 

0.00 

0.00 

[60] 

115 

3CgHn 

+  H02 

= 

C8H9  +  CH20  +  OH 

2.5000E+09 

0.00 

0.00 

[60] 

116 

3CgHn 

= 

c8h8  +  CH3 

2.0000E+13 

0.00 

113.10 

[60] 

117 

1C9Hh 

= 

c8h8  +  ch3 

1.0000E+14 

0.00 

145.81 

[96] 

118 

'CsHu 

+  O 

= 

c7h6o  +  c2h5 

1.6000E+10 

0.00 

0.00 

[96] 

119 

'CsHu 

+  OH 

= 

C7H60  +  C2H6 

1.6000E+10 

0.00 

0.00 

[96] 

120 

2C9Hh 

+  OH 

= 

C7H8  +  CH3CHO 

2.0000E+10 

0.00 

16.73 

[96] 

121 

2C9Hh 

= 

C6H5  +  C3H6 

1.0000E+14 

0.00 

145.60 

[96] 

122 

3CgHn 

= 

C7H7  +  C2H4 

1.2300E+13 

-0.10 

126.37 

Present 

Work 

123 

CgHi2 

+  H 

= 

=  1CgHn  +  H2 

1.0100E+01 

2.90 

8.07 

[96] 

124 

CgHi2 

+  O 

= 

=  1C9Hh  +  OH 

4.7700E+01 

2.70 

4.62 

Adj.[1 04] 
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No 

Reaction 

A 

n 

Ea 

Ref 

125 

C9H12 

+  OH 

=  1CgHi1  +  H20 

7.0800E+03 

1.90 

-4.84 

[96] 

126 

C9H12 

+  H 

=  2C9Hh  +  H2 

1.0100E+01 

2.90 

12.26 

[96] 

127 

C9H12 

+  0 

=  2C9Hh  +  OH 

4.7700E+01 

2.70 

8.81 

Adj.[104] 

128 

C9H12 

+  OH 

=  2C9Hn  +  H20 

4.7865E+03 

1.40 

3.56 

[60] 

129 

C9H12 

+  H 

=  3C9Hn  +  H2 

1.3300E+03 

2.50 

28.26 

[96] 

130 

C9H12 

+  0 

=  3C9Hh  +  OH 

1.9300E+02 

2.70 

15.54 

[96] 

131 

C9H12 

+  OH 

=  3C9Hh  +  H20 

3.1600E+04 

1.80 

3.90 

Adj.[212] 

132 

C9H12 

+  H02 

=  1C9Hh  +  H2O2 

9.6400E+00 

2.60 

54.01 

[96] 

133 

C9H12 

+  ch3 

=  1C9Hh  +  CH4 

1 .5000E-03 

3.46 

22.90 

Adj.[104] 

134 

C9H12 

+  H02 

=  2C9Hh  +  H2O2 

9.6400E+00 

2.60 

58.20 

Adj.[104] 

135 

C9H12 

+  H02 

=  3CgHn  +  H2O2 

4.7600E+01 

2.50 

69.01 

Adj.[104] 

136 

C9H12 

+  ch3 

=  3C9Hh  +  CH4 

9.0300E-04 

3.65 

29.90 

Adj.[104] 

137 

C9H12 

+  ch3 

=  2C9Hh  +  CH4 

1 .5055E-03 

3.46 

22.90 

Adj.[104] 

138 

C9H12 

+  C2H5 

=  ^sHn  +  C2H6 

1.2100E-03 

3.46 

31.26 

Adj.[1 04] 

139 

C9H12 

+  C2H5 

=  2C9Hn  +  C2H6 

1.2100E-03 

3.46 

31.26 

Adj.[104] 

140 

C9H12 

+  C2H5 

=  3CgHn  +  C2H6 

9.0400E-04 

3.65 

38.24 

Adj.[104] 

141 

C9H12 

+  C2H3 

=  ^gHu  +  C2H4 

1.0000E+00 

3.10 

36.94 

Adj.[104] 

142 

C9H12 

+  c2h3 

=  2C9Hh  +  C2H4 

1.0000E+00 

3.10 

36.94 

Adj.[104] 

143 

C9H12 

+  c2h3 

=  3C9Hn  +  C2H4 

6.0000E-01 

3.30 

43.94 

Adj.[1 04] 

144 

C9H12 

+  c3h5(A) 

=  ^sHn  +  C3H6 

7.9400E+08 

0.00 

67.82 

[96] 

145 

C9H12 

+  c3h5(A) 

=  2C9Hh  +  C3H6 

7.9400E+08 

0.00 

67.78 

[96] 

146 

C9H12 

+  c3h5(A) 

=  3C9Hh  +  C3H6 

7.9400E+08 

0.00 

85.77 

[96] 

147 

C9H12 

+  c6h5 

=  ^gHn  +  C6H6 

7.9400E+08 

0.00 

67.78 

[96] 

148 

C9H12 

+  C6H5 

=  2C9Hh  +  C6H6 

7.9400E+08 

0.00 

67.68 

[96] 

149 

C9H12 

+  c6h5 

=  3C9Hh  +  C6H6 

7.9400E+08 

0.00 

85.77 

[96] 

150 

C9H12 

+  C7H7 

=  ^gHu  +  C7H8 

7.9400E+08 

0.00 

67.78 

[96] 

151 

C9H12 

+  C7H7 

=  2C9Hh  +  C7H8 

7.9400E+08 

0.00 

67.68 

[96] 

152 

C9H12 

+  C7H7 

-  3C9Hn  +  C7H8 

7.9400E+08 

0.00 

85.84 

[96] 

153 

C9H12 

+  02 

=  ^gHii  +  HO2 

4.0000E+10 

0.00 

149.78 

[96] 

154 

C9H12 

+  02 

=  2C9Hii  +  H02 

4.0000E+10 

0.00 

207.10 

[96] 

155 

C9H12 

+  02 

=  3C9Hh  +  H02 

4.0000E+10 

0.00 

211.71 

[96] 

156 

C9H12 

=  C6H5  +  C3H7(N) 

1.0000E+16 

0.00 

305.43 

Adj.[13] 

157 

C9H12 

=  C7H7  +  C2H5 

1.0000E+16 

0.00 

305.43 

Adj.[13] 

158 

C9H12 

=  CsHg  +  CH3 

8.0000E+16 

0.00 

353.13 

Adj.[13] 

159 

C9H12 

+  H 

=  C3H7(N)+  c6h6 

5.7800E+10 

0.00 

33.84 

Adj.[80] 
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Reaction 

A 

n 

Ea 

Ref 

160 

C9H12  =  3C9Hh  +  H 

8.0000E+15 

0.00 

419.27 

[96] 

161 

C9H12  =  2C9Hh  +  H 

8.0000E+15 

0.00 

414.00 

[96] 

162 

C9H12  =  ^gHii  +  H 

8.0000E+15 

0.00 

360.70 

[96] 

163 

C9H12  +  C3H5(T)  =  ^gHn  +  C3H6 

7.9400E+08 

0.00 

67.78 

[96] 

164 

C9H12  +  C3H5(S)  =  ^sHn  +  C3H6 

7.9400E+08 

0.00 

67.78 

[96] 

165 

C9H12  +  C3H5(T)  =  2C9Hh  +  C3H6 

7.9400E+08 

0.00 

67.78 

[96] 

166 

C9H12  +  C3H5(S)  =  2C9Hh  +  C3H6 

7.9400E+08 

0.00 

67.78 

[96] 

167 

C9H12  +  C3H5(T)  =  3C9Hh  +  C3H6 

7.9400E+08 

0.00 

85.77 

[96] 

168 

C9H12  +  C3H5(S)  =  3C9Hn  +  C3H6 

7.9400E+08 

0.00 

85.77 

[96] 

Enhanced  Collision  Efficiencies 

Unless  otherwise  specified  all  collision  efficiencies  are  set  equal  to  1. 


